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The . program  provides  a  uniform  method  for  calculating  wave  runup  and  wind  set¬ 
up  at  the  edge  of  lakes  and  reservoirs.  Elevations  of  wave  runup  can  be 
computed  for  embankment  slopes  that  are  either  smooth,  turfed,  or  riprapped. 
References:  Refer  to  Appendix  A. 
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The  program  is  written  in  G635  time-share  series,  FORTRAN  IV,  and  is  part  of 
a  Conversationally  Oriented  Real-Time  Program-Generating  System  (CORPS). 


EQUIPMENT  OCT  AILS 


The  basic  program  was  developed  on  CDC  7600,  Southwestern  Division,  Dallas,  TX. 
The  revised  computer  version  in  this  abstract  is  now  operational  on  the  WES 
G635,  Vicksburg,  MS;  HIS  66/80,  Macon,  GA;  and  Boeing  CDC,  Seattle,  WA. 


All  inputs  for  program  H7780  are  either  read  from  a  permanent  data  file  or 
are  cued  and  read  upon  entry  to  the  input  subroutine  H7780I.  All  specific 
input/output  requirements  for  program  H7780  are  given  in  PART  II:  COMPUTER 
FUNCTIONAL  DESCRIPTION  of  this  abstract. 


The  basic  program  was  developed  by  B.  R.  Bodine,  primarily  bpsed  on  the  infor¬ 
mation  presented  in  ETL  1110-2-221,  "Wave  Runup  and  Wind  Setup  on  Reservoir 
Embankments"  by  Bruce  L.  McCartney.  The  coding  has  been  revised  at  WES  in 
order  to  accommodate  the  CORPS  time-share  feature^.  i 
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PROGRAM  SPECIFICATIONS: 

Language :  ANSI  FORTRAN  (FORTRAN  IV) 

Solution  Requirements:  The  run  command 

RUN  WESLIB/CORPS/H7780,  R 
and  all  necessary  inputs. 

Method  of  Analysis:  Wave  runup  and  wind  setup  are  solved  by  direct 
solution  of  algebraic  equations  except  for  the  transitional  depth 
wave  length  which  is  solved  via  an  interactive  technique. 

Core  Requirements  G635:  J  5  K  words 

External  Storage:  Type  of  storage  -  disk.  Each  user  may  have  an 
existing  input  data  file  with  a  maximum  of  12  K  G635  ASC  II  words. 

If  the  file  does  not  exist,  a  user  named  file  is  created.  If  user 
decides  to  save  tabulated  output  for  graphics  and/or  other  use,  each 
output  file  may  have  a  maximum  of  12  K  G635  ASC  II  words.  This  file 
may  already  exist  or  be  created  by  the  user  during  a  run  of  H7780. 
Restrictions :  Refer  to  pages  18-22 

General  Equations :  Refer  to  equations  10-18  as  presented  on 
pages  lh-17. 
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REF:  ER  1110-1-10  -  ENGINEERING  AND  DESIGN  -  Engineering  Computer 
Program  Library  Standards  and  Documentation,  Appendix  B 

PART  I:  ENGINEERING  DESCRIPTION 

1.  PROGRAM  NUMBER:  723-X6-M007A 

2.  TITLE :  H7780  -  Wave  Runup  and  Wind  Setup  -  Computational  Model 

3.  REVISION  LOG:  Program  coding  revised  June  j8  to  accommodate  the 
CORPS  time-share  features. 

4.  PURPOSE  OF  PROGRAM:  To  provide  a  uniform  method  for  calculating 
wave  runup  and  wind  setup  at  the  edge  of  lakes  and  reservoirs.  Eleva¬ 
tions  of  wave  runup  can  be  computed  for  embankment  slopes  that  are 
either  smooth,  turfed,  or  riprapped. 

References:  Refer  to  Appendix  A. 

5.  STEP  SOLUTION: 

a.  Introduction: 

The  computer  program  outlined  in  this  abstract  has  been  developed 
specifically  for  estimating  wave  runup  and  wind  setup  at  the  edge  of 
Lakes  or  reservoirs.  This  program  may  also  be  used  to  estimate  the 
elevations  of  wave  runup  at  the  edge  of  semienclosed  water  bodies  such 
as  bays  and  estuaries.  However,  the  method  used  in  the  present  solution 
scheme  for  estimating  wind  setup  is  not  generally  appropriate  for  semi¬ 
enclosed  water  bodies.  The  principles  involved  and  estimating  pro¬ 
cedures  used  in  developing  the  program  are  presented  together  with  the 
limitations  of  the  computational  system.  Example  problems  involving  the 
determination  of  wave  runup  and  wind  setup  in  more  or  less  typical  lakes 
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or  reservoirs  are  presented  to  demonstrate  the  utility  of  the  program 
and  procedures  necessary  for  application. 

This  abstract  is  not  intended  to  stand  on  its  own  but  rather  serve 
as  a  supplement  to  ETL  1110-2-221  (McCartney,  1976)  references  in 
Appendix  A.  Inasmuch  as  the  referenced  ETL  is  self-contained  in  respect 
to  estimating  characteristics  of  wind-generated  waves  for  deepwater 
conditions,  guidance  for  estimating  waves  for  transitional  water  and 
shallow  water  conditions  is  referenced  to  the  Shore  Protection  Manual 
(1977).  In  the  determination  of  freeboard  allowances  for  high  dams, 
water  depths  are  usually  deep  and  thus  it  is  usually  necessary  to  esti¬ 
mate  wave  runup  on  waves  generated  in  deep  water  or  nearly  deep  water. 
However,  in  some  instances,  the  design  of  some  slope  protection  mea¬ 
sures,  such  as  highway  and  railroad  embankments  and  flood  levees,  the 
water  depths  over  the  wave  generating  area  may  be  small  and  require  that 
the  wave  characteristics  be  based  on  transitional  water  or  shallow  water 
conditions.  Consequently,  the  present  computational  scheme  includes 
wave  prediction  for  both  deep,  transitional,  and  shallow  water. 

The  program  has  been  formulated  in  such  a  manner  that  the  prepara¬ 
tion  of  basic  data  is  minimized.  In  addition,  provisions  have  been  made 
to  allow  determinations  of  wave  runup  and  wind  setup  for  one  or  more 
sites  with  a  single  computer  run. 

b.  Input :  The  inputs  are  read  from  existing  user  input  data  file 
or  are  entered  at  run  time  via  the  user's  input  terminal  device.  The 


2 


H7780 


inputs  are  described  in  PART  II:  COMPUTER  FUNCTIONAL  DESCRIPTION, 

Part  B,  INPUT  DATA  DESCRIPTION. 

c .  Mathematical  Formulation: 

1.  Design  Wind.  In  the  present  program  two  options  are  provided  for 
describing  the  design  wind.  The  first  of  these  is  simply  the  specifi¬ 
cation  of  the  design  wind  for  a  particular  location  based  on  actual  wind 
records  representative  of  the  area.  For  this  option  the  overland  wind 
speed,  U  ,  in  miles  per  hour  is  a  basic  input  to  the  program.  The 
overland  wind  speed  is  converted  to  the  wind  speed  over  water  internal 
to  the  program  based  on  previously  adopted  criteria  presented  in 
ETL  1110-2-221.  These  criteria  relate  the  ratio,  R  ,  wind  speed  over 
water  and  land  areas  to  the  effective  fetch,  F  ,  for  wave  generation 
as  follows. 

F  (miles)  0.5  1  2  3  ^  5  (or  over) 

R^  =  ^oyerwater  1>0g  1>13  1>2i  1>26  1>2g  1>30 

overland 

For  convenience  in  programing  the  relationship,  the  ratio,  R^  ,  is 
presumed  to  approach  unity  as  the  fetch  approaches  zero. 

In  the  absence  of  sufficient  wind  observations  for  a  particular 
site,  the  design  wind  may  be  determined  from  the  averaged  maximum 
regional  wind  data  provided  in  ETL  1110-2-221.  For  this  case,  the  maxi¬ 
mum  1-minute  and  60-minute  averaged  winds  for  a  specific  location  are 
determined  from  the  appropriate  figures  (Figures  2  through  9)  in  the  ETL. 
These  overland  wind  speeds  are  a  basic  input  to  the  program.  The 
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adjustment  to  overwater  wind  speeds  are  determined  internal  to  the 
program.  Moreover,  the  design  wind  speed  is  calculated  by  a  program 
subroutine.  The  techniques  and  procedures  used  for  calculating  the 
design  wind  speed  is  discussed  in  Appendix  B. 

2.  Effective  Fetch  for  Wave  Generation.  In  the  present  computational 
scheme  the  effective  fetch  distance  for  wave  generation  is  determined 
in  accordance  to  the  method  developed  Oj  Saville  et  al.  (1962)  and  out¬ 
lined  in  ETL  1110-2-221.  The  procedure  for  determining  the  effective 
fetch  distance,  F  ,  by  this  method  is  illustrated  in  Figure  1.  It  con¬ 
sists  of  superimposing  radials  on  a  suitable  map  which  depicts  the  water 
body.  The  radials  are  constructed  in  a  manner  so  that  they  emanate  from 
the  shoreline  site  where  wave  information  is  needed  and  extend  across 
the  water  area  until  they  intersect  the  shoreline.  Even  though  some  of 
the  radial  lengths  may  be  zero,  the  central  radial  taken  along  the  long 
axis  of  the  water  body  will  usually  result  in  the  longest  effective  wave 
fetch.  The  interval  between  any  two  adjacent  radials  is  taken  as 
6  degrees  (angle  denoted  by  the  symbol,  a  ,  on  Figure  l)  and  the 
radials  are  limited  by  an  angle  of  U5  degrees  on  either  side  of  the 
wind  direction.  As  a  result  of  the  bounds  indicated,  15  radials  are 
used  in  the  determination  of  the  effective  fetch. 

In  the  application  of  the  effective  fetch  method  in  practice,  ques¬ 
tion  sometimes  arises  as  to  whether  the  radials  have  been  constructed 
properly  for  maximum  wave  generation.  On  this  basis,  provisions  have 
been  made  in  the  present  computational  scheme  to  input  more  than  15 
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radials,  if  desired.  A  special  routine  is  included  in  the  program  to 
determine  the  effective  fetch  with  the  greatest  length  (referred  to 
hereafter  as  the  critical  effective  wave  fetch)  from  a  particular  set  of 
15  sequential  radials.  It  is  to  be  noted,  however,  that  a  minimum  of 
15  radials  are  required  as  a  basic  input  to  the  program. 

3.  Effective  Fetch  for  Wind  Setup.  Fetch  distance  for  use  in  estimat¬ 
ing  wind  setup  can  usually  be  taken  considerably  longer  than  fetch  dis¬ 
tance  used  for  wave  generation.  This  is  because  wind  setup  effects  may 
be  transferred,  to  some  extent,  around  substantial  bends  in  the  water 
body,  thus  warranting  the  assumption  of  a  longer  fetch.  Guidance  pro¬ 
vided  in  ETL  1110-2-221  indicates  that  the  wind  setup  fetch,  ,  can 
be  estimated  by  using  twice  the  effective  wave  fetch.  For  most  problems 
encountered  in  practice,  this  is  a  sufficient  estimate  of  the  wind  setup 
fetch.  However,  for  some  problems  such  an  estimate  may  not  always  be 
realistic.  Consider  the  example  that  all  radial  lengths  are  approxi¬ 
mately  equal  for  a  given  wave  generation  area.  This  would  imply,  as  can 
be  rhown,  that  the  effective  wave  fetch  is  approximately  equal  to  88 
percent  of  the  central  radial  length.  Because  of  this  particular  con¬ 
dition,  the  wind  setup  fetch  would  exceed  the  central  radial  length  by 
about  76  percent  when  using  twice  the  effective  wave  fetch.  As  a  conse¬ 
quence,  this  would  infer  that  a  portion  of  the  wind  setup  fetch  would  be 
over  water  and  another  portion  of  the  fetch  would  be  over  land.  In  the 
present  computational  scheme,  this  possibility  is  avoided  by  restricting, 
the  wind  setup  fetch  to  a  distance  less  than  or  equal  to  the  longest 
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radial  length.  Specifically,  the  wind  setup  fetch  is  equal  to  twice  the 
effective  wave  fetch  if  this  distance  is  less  than  the  longest  radial 
and  equal  to  the  longest  radial  if  2  F  is  longer  than  the  longest 
radial . 

1,1  Wave  Character i sties.  The  most  fundamental  description  of  waves  in¬ 
duced  by  wind  is  their  length  L  ,  height  I[  ,  period  T  ,  and  the  wave 
depth  d  (see  Figure  2).  Inasmuch  as  wind  waves  are  merely  a  more  or 
less  rhythmic  change  in  the  elevation  of  the  water  surface  with  essen¬ 
tially  no  net  transport  of  the  water,  the  wave  form  propagates  in  the 
direction  of  the  Wind  by  a  velocity  equal  to  L/T  ,  This  velocity  is 
referred  to  as  the  phase  velocity  or  wave  celerity. 

Surface  waves  induced  by  wind  are  classified  in  accordance  to  their 
lengths  and  the  depths  of  water  over  which  they  travel,  or  specifically 
by  ratio  d/L  .  A  wave  is  said  to  be  a  "deep  water  wave"  if  d/I,  is 
greater  than  j./ 2 ,  a  transitional  wave"  if  d/L  lies  in  a  range  between 
1/25  and  1/2,  and  a  "shallow  water  wave"  if  d/L  is  less  than  1/25. 

The  characteristics  of  deepwater  waves  are  unaffected  by  the  depth  over 
which  they  propagate.  On  the  other  hand,  waves  in  transitional  and 
shallow  water  are  affected  by  the  depth  in  which  they  run,  resulting 
primarily  in  the  modification  of  the  wave  form.  It  is  often  said  that 
transitional  and  shallow  water  waves  feel  bottom. 

Observing  waves  under  the  influence  of  wind  reveal  that  there  is 
little  regularity  in  wave  form  and  the  elevations,  of  the  rise  and  fall 
of  the  water  surface.  When  attention  is  focused  on  say  1,000 
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consecutive  waves,  referred  to  as  a  wave  train,  passing  a  fixed  location 
it  will  be  seen  that  there  is  considerable  irregularity  in  wave  form  and 
wave  height.  Once  the  waves  leave  the  generating  area,  the  waves  con¬ 
tinue  to  propagate  in  the  direction  the  wind  was  blowing,  but  they  begin 
to  transform  into  waves  of  more  regular  shape  and  more  uniform  height. 

In  the  solution  of  engineering  problems  for  reservoirs  and  lakes,  wave 
analysis  is  almost  always  concerned  with  waves  under  the  influence  of 
wind.  Thus  in  dealing  with  such  problems  it  is  necessary  from  a  theo¬ 
retical  point  of  view  to  represent  the  waves  by  an  idealized  description 
of  the  water  surface.  To  facilitate  this  description  a  wave  referred  to 
as  the  "significant  wave"  has  been  postulated  in  which  this  wave  has 
been  almost  universally  adopted  for  the  solution  of  wave  problems.  The 
significant  wave  height  is  defined  as  the  average  height  of  the  one- 
third  highest  waves  and  written  as  Hs  or  •  The  wave  period 

associated  with  such  a  wave  is  called  the  significant  wave  period,  T^  . 

Occasionally,  in  the  planning  and  design  of  shore  structures  it  is 
advantageous  to  consider  other  waves  in  the  wave  spectrum  for  particular 
design  analysis.  This  stems  from  the  fact  that  in  some  instances  it  may 
be  necessary  or  appropriate  to  provide  a  higher  or  lesser  degree  of  pro¬ 
tection  than  that  afforded  by  the  significant  wave.  Altering  the  degree 
of  protection  may  be  accomplished  by  selecting  a  specific  wave  height 
greater  or  smaller  than  the  significant  wave  height.  Other  wave  heights 
may  be  obtained  and  based  on  a  knowledge  of  the  significant  wave  height 
and  utilizing  the  Rayleigh  distribution  function.  The  wave  height 
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distribution,  accordingly,  is  given  by 

Hp  .  Ks(^  ) 


where  is  the  wave  height  associated  with  a  particular  probability 

of  exceedence,  p  ,  and  In  signifies  the  natural  logarithm.  For 
example,  a  1  percent  wave  height  (p  =  0.01)  i.e.,  the  elevation  exceeded 
by  1  percent  of  the  wave  heights,  according  to  equation  (l)  yields 


K0.01  =  Hs 


In  (1/0.01) 


1/2 

=  1.517  H 

s 


Thus  the  1  percent  wave  is  about  1.5  times  larger  than  the  significant 
wave.  The  present  solution  scheme  allows  selection  of  a  wave  other  than 
the  significant  wave  for  design  purposes. 

5.  Wave  Forecasting.  A  great  deal  of  research  effort  has  been  devoted 
to  the  prediction  of  waves  under  the  influence  of  the  wind  that  gene¬ 
rated  them.  As  a  result  numerous  wave  prediction  schemes  have  become 
available.  Some  of  these  schemes,  are  complex  while  others^ are  rather 
simple.  A  relatively  simple  wave  prediction  scheme  which  has  been  used 
rather  extensively  in  practice  is  called  the  Sverdrup-Munk-Bretschneider 
(SMB)  method.  This  method,  as  revised  by  Bretschneider  in  1971 »  is  dis¬ 
cussed  in  the  Shore  Protection  Manual  (SPM,  1977)  and  is  used  for  the 
method  of  wave  forecasting  in  connection  with  the  presexit  solution 
scheme.  The  method  involves  estimating  the  wave  characteristics  (height 
and.  period)  and  the  time  period,  called  the  duration,  necessary  to 
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generate  the  wave  by  semiempirical — semitheoretical  formulae  or  by  avail¬ 
able  wave  forecasting  curves  (see  SPM,  1977).  The  formulae,  used  in 
this  development  is  presented  in  the  following  paragraphs  for  deep, 
transitional,  and  shallow  water. 

a.  Deepwater  waves.  The  wave  forecasting  equations  for  deep  water 

are , 

(2) 


Si 

u2 


0.283  tan  h 


°-0125(P) 


0.1*2 


where 

,  ,  {x} 

exp  (x>  =  e 

In  =  log 

e 

K  =  6.5882 
A  =  0.0161 
B  =  0.3692 
C  =  2.202U 

and 

D  =  0.8798 

The  notations  used  are  as  follow:  g  is  the  acceleration  of  gravity; 

H  is  the  wave  height;  U  is  the  wind  speed;  F  is  the  effective  wave 
fetch  distance;  t  is  the  duration;  T  is  the  wave  period;  and  tanh 
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is  the  hyperbolic  tangent.  Equations  2,  3,  and  1  are  nondimensional ; 
and  consequently,  any  system  of  units  may  be  used.  Units  consisting  of 
feet-seconds  are  used  in  the  present  computational  scheme. 

Subsequent  to  the  determination  of  the  wave  period,  the  wave  length 
in  deep  water  Lq  may  be  calculated  by 

2 

Lq  =  =  5-12  T2  (ft)  (5) 


This  relation  shows  that  the  length  of  the  wave  in  deep  water  is  unaf¬ 
fected  by  the  depth  in  which  the  wave  travels  (i.e.,  where  d/LQ  >  0.5). 

b.  Transitional  and  shallow  water  waves.  Wave  generation  and  wave 
characteristics  are  affected  when  waves  are  traveling  in  transitional 
and  shallow  water  depths.  For  a  given  set  of  wind  and  fetch  conditions, 
wave  heights  will  be  smaller  and  wave  periods  shorter  in  comparison  to 
those  in  deep  water. 

The  numerical  technique  presented  here  is  based  on  successive 
approximations  in  which  wave  energy  is  added  due  to  wind  stress  and 
subtracted  due  to  bottom  friction  and  percolation. 

The  equations  applicable  for  transitional  and  shallow  water  condi¬ 
tions  which  incorporate  bottom  friction  are  as  follows: 


=  0.283  tanh 

U 


0.530 


\0.75 


0.0125 


f£ 


,0.U2 


tanh 


[  tanh 


0.5301 


gd 


,0.75 


(6) 
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and 


_ST 


2nU 


=  1.20  tanh 


0. 833(^) 

\U  / 


0.375 


tanh 


(  o.«n(£ 

j  0.25 

/  tanh 

0.833( 

k\°-37? 

V 

OJ2/ 

(7) 


in  which  in  deep  water  reduce  to  equations  (2)  and  (3).  Equations  (6) 
and  (7)  are  considered  valid  provided  that  wind  blows  over  the  water 
with  sufficient  duration  to  fully  develop  the  waves. 

In  shallow  water  the  wave  length,  L,  is 


L  =  T  */gd  (8) 

in  which  the  term  i/g(l  is  the  wave  celerity  in  shallow  water. 

For  transitional  depths  the  wave  length,  L,  is  given  by 

1  ■  Sr  '*>{¥)  <»> 

Evaluation  of  the  wave  length  from  equation  (9)  involves  some  diffi¬ 
culty  since  the  unknown,  L  ,  appears  on  both  sides  of  the  equation. 

In  the  computer  program  developed,  a  special  iterative  and  convergent 
scheme  is  used  in  the  determination  of  L  . 

6.  Wave  Runup.  The  vertical  height  above  the  Stillwater  level  (SWL) 
that  a  wave  will  run  up  the  face  of  a  structure  {see  Figure  2)  depends 
on  several  factors.  These  factors  are  identified  as  structure  shape 
and  roughness,  water  depth  at  structure,  bottom  slope  in  front  of  the 
structure,  and  the  characteristics  of  the  waves  impinging  on  the 
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structure.  Because  of  the  large  number  of  variables  involved  and  possi¬ 
ble  combinations  of  geometric  shapes  and  wave  conditions,  an  accurate 
description  of  wave  runup  phenomena  is  not  available.  However,  a  great 
deal  of  useful  guidance  has  become  available  for  estimating  wave  runup 
as  a  result  of  numerous  laboratory  studies  conducted  by  using  mechanical 
generated  waves  and  some  field  studies.  Seville  (1955>  1956)  Savage 
(1959),  Saville  et  al.  (1962),  Battjes  (197*+),  and  Ahrens  and  McCartney 
(1975),  among  others,  present  the  results  of  such  studies. 

A  large  number  of  laboratory  experiments  were  conducted  for  runup 
on  smooth  plane  slopes  in  1956  (Saville,  1956).  The  results  of  these 
experiments  were  reproduced  in  graphic  form  and  are  presented  in  the 
recent  Shore  Protection  Manual  (1977)-  In  1965,  Franzius,  a  German 
Engineer,  transformed  these  graphical  runup  results  into  formulae  which 
are  convenient  for  computer  applications.  The  relations  obtained,  as 
given  by  Battjes  (197*+),  are  as  follows: 


R  =  H  sin  9  (5-95  tan  9  +  1,5) 


(m) 


(10) 


where  the  exponent  a  is  given  by 

.58  -  2.35  tan  0)  +  0.092  cot  9  -  0.26 


A“-: 


(11) 


in  which  R  is  the  wave  runup,  ds  is  the  depth  of  the  water  at  the 
toe  of  the  embankment,  and  9  is  the  angle  between  the  horizontal  and 
the  embankment  slope.  Equation  (10)  is  valid  for 
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1/6  <  tan  e  <  1/2.25  and  £  O.U75 

s 

For  a  1  on  1.5  slope  the  following  equation  is  given 


(12) 


R  =  2.3  H 


-  0.18) 


(13) 


For  slopes  flatter  than  a  1  on  6  slope,  a  formula  given  Battjes 
(1975)  is  used  to  estimate  the  wave  runup.  This  equation  is  Riven  by 


R  =  0.1  T  •'gH  tan  8 


1 1 
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Franz ius  also  found  from  experiments  that  the  runup  on  smooth  slopes 
could  be  modified  to  account  for  turfed  slopes  by  the  following  approxi¬ 
mate  relation 


=  (0.85  to  0.90)  R  (15) 

I  s 

in  which  Rm  is  the  runup  on  a  turfed  slope  and  R  is  the  runup  on  a 
1  s 

smooth  slope.  In  other  words,  the  runup  on  a  turfed  slope  is  about  85 
to  90  percent  of  the  runup  on  a  smooth  slope.  An  approximate  relation 
similar  to  equation  (15)  may  also  be  Riven  for  impermeable  stepped 
slopes  in  which  the  vertical  and  horizontal  dimensions  of  the  steps 
are  relatively  small  compared  to  the  wave  height.  Such  a  relation  ir 


R_  =  (0.7  to  0.8)  R 
1  £» 


(16) 
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The  factors  0.7  to  0.8  were  based  on  data  presented  in  SPM  (1977)  for 
stepped  slopes  and  private  communications  with  John  Ahrens,  Oceanogra¬ 
pher,  Coastal  Structures  Branch,  Coastal  Engineering  Research  Center. 

For  solution  of  other  engineering  problems,  it  is  reasoned  that 
runup  for  smooth  slopes  can  be  modified  in  the  same  manner  for  other 
slope  surface  conditions  provided  that  the  slope  is  uniform  and  the 
embankment  surface  is  relatively  smooth.  Such  a  modification  would  not, 
however,  be  applicable  to  extremely  rough  slope  surfaces. 

In  the  case  of  riprapped  embankment  slopes,  the  relation  given  by 
McCartney  (1976)  provides  a  reasonable  estimate  of  the  runup.  This 
relation  is  given  by 


0.U  +  (H/Lo)1//2  cot  6 


Equation  (17)  is  considered  valid  for  2  <  cot  6  <  5  . 

The  equations  given  herein  for  estimating  runup  will  normally  be 
sufficient  for  resolving  most  problems  encountered  in  reservoirs  and 
lakes . 

7.  Wind  Setup.  The  action  of  wind  blowing  over  water  surfaces  of  lakes 
and  reservoirs  is  not  only  responsible  for  generating  surface  waves  but 
causes  the  water  surface  to  tilt  from  the  windward  side  to  the  leeward 
side  of  the  basin.  The  tilt  is  a  result  of  the  wind  inducing  a  current 
in  the  upper  layers  of  the  water  in  the  direction  of  the  wind  and  thus 
causing  the  water  to  pile  up  at  the  leeward  shore.  A  return  current. 
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smaller  in  magnitude  than  the  wind  induced  surface  current,  is 
established  along  the  bottom  in  the  opposite  direction  of  the  wind. 

At  the  upwind  side  of  the  basin  the  water  level  is  depressed  and  at  the 
downwind  side  the  water  level  is  raised  as  a  result  of  the  unequal  sur¬ 
face  and  bottom  current.  The  rise  in  water  level  at  the  downwind  side 
of  the  water  body  is  referred  to  as  wind  setup.  The  prediction  of  wind 
setup  is  complex  due  to  the  number  of  mechanisms  and  processes  involved 
and  the  complications  that  arise  in  simulating  the  basin  geometry.  Al¬ 
though  wind  setup  can  be  estimated  with  a  reasonable  degree  of  accuracy 
by  sophisticated  mathematical  models  presently  available,  there  is 
usually  little  justification  to  do  so  for  relatively  small  and  deep 
lakes  and  reservoirs  due  to  the  substantial  effort  and  expense  in  apply¬ 
ing  such  schemes  and  the  fact  that  wind  setup  is  rather  small  in  coin- 
par  ison  to  the  wave  setup. 

A  relation  which  in  general  gives  a  reasonably  good  estimate  of 
the  wind  setup  (see  McCartney,  1976)  is 


1UU0  D 


in  which  S  is  the  wind  setup  in  feet;  U  is  the  wind  speed  in  miles 
per  hour;  F  is  the  wind  setup  fetch  in  miles  and  D  is  the  average 
depth  over  the  wind  setup  fetch  in  feet. 

Equation  (l8)  will  generally  provide  conservative  wind  setup  esti¬ 
mates  at  all  locations  in  a  relatively  deepwater  body  except  at  those 
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locations  where  the  basin  geometry  converges  to  the  shore  site  where 
interest  is  centered.  For  sites  located  in  a  highly  convergent  zone  it 
may  be  justified  to  increase  the  wind  setup  as  much  as  50  percent  to 
account  for  an  additional  pileup  of  water. 

d.  Output :  The  output  consists  of  three  parts;  basic  input,  com¬ 


putational  results,  and  user  option  output.  The  basic  input  and  compu¬ 
tational  results  are  printed  for  each  problem.  The  user  option  output 
(Wind  Data  and  Deepwater  Wave  Characteristics  and  Wave  Heights  and  Wave 
Runup  for  Waves  Other  Than  the  Significant  Wave)  is  or  is  not  printed 
according  to  user  specification  in  the  basic  input.  The  tabulated  out¬ 
put  (user  option  output),  if  desired,  may  be  stored  in  a  user  existing 
or  run-time  created  output  data  file. 

6.  ACCURACY :  Governed  by  accuracy  of  input  data. 

7.  REMARKS: 

a.  General ■  Inasmuch  as  the  computational  model  described  herein 
will  usually  provide  a  reasonable  estimate  of  wave  runup  and  wind  setup 
for  most  problems  encountered  in  practice,  there  are  some  practical 
problems  in  which  the  model  cannot  be  expected  to  yield  meaningful 
results.  Because  the  model  is  not  applicable  to  all  problems,  care 
must  be  exercised  in  determining  the  practical  limits  in  which  the  model 
can  be  applied.  Determination  of  whether  the  model  is  applicable  or 
not  applicable  to  a  specific  problem  area  requires  a  careful  review  of 
the  physical  characteristics  of  the  water  body  and  embankment  coupled 
with  engineering  judgment. 
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The  model  is  limited  to  the  determination  of  wave  runup  on  uniform 
smooth,  relatively  smooth,  and  riprap  slopes  with  waves  approaching  from 
a  direction  normal  (90  degrees)  to  the  slope.  It  is  also  applicable  to 
stepped  slopes  in  which  the  vertical  and  horizontal  dimensions  of  the 
steps  are  small  in  comparison  to  the  wave  height.  The  model  is  not  ap¬ 
plicable  to  other  nonuniform  slope  configurations,  vertical  walls,  and 
for  such  slope  covers  as  permeable  rubble  mound.  Although  the  model 
does  not  apply  for  the  cases  mentioned,  wave  runup  estimates  can  usually 
be  made  by  the  methods  and  procedures  presented  in  SPM  (1977). 

The  height  and  period  of  waves  calculated  in  the  model  are  deter¬ 
mined  from  empirical  wave  forecast  formulae  in  which  it~is  assumed  that 
wave  propagation  over  the  fetch  occurs  with  a  uniform  depth.  However, 
for  real  water  bodies,  depths  can  vary  considerably  over  the  wave  fetch. 
In  many  cases  when  waves  propagate  over  the  fetch,  the  wave  in  one  por¬ 
tion  of  the  fetch  may  be  a  deepwater  wave  while  at  another  a  shallow 
water  wave.  However,  a  reasonable  estimate  of  the  height  and  period  of 
waves  can  be  obtained  by  using  the  average  depth  over  the  fetch  provided 
that  the  wave  does  not  break  as  it  travels  over  the  fetch.  Theoretic¬ 
ally,  a  wave  height,  H  ,  cannot  exceed  0.78d,  where  d  is  the  local 
water  depth.  Waves  breaking  offshore  from  a  structure  result  in  a  lower 
wave  height  at  the  structure.  Thus,  the  reduced  wave  height  should  be 
used  in  determining  the  wave  runup  at  the  structure.  No  provisions 
are  made  in  the  present  model  to  account  for  breaking  waves  and,  there¬ 
fore,  it  is  necessary  to  account  for  this  effect  separately. 
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There  are  some  instances  when  waves  break  along  the  path  of  wave 
propagation  at  a  considerable  distance  from  the  site  where  interest  is 
centered.  For  example,  a  roadway  may  transverse  the  water  body  at  about 
right  angles  to  the  wind  direction  at  about  mdiway  along  the  wave  fetch. 
Under  severe  flooding  conditions,  the  roadway  may  cause  the  waves  to 
break  as  they  pass  over  the  top  of  the  road.  For  such  cases,  an  esti¬ 
mate  of  the  design  wave  reaching  the  project  site  may  be  obtained  by 
first  obtaining  the  breaking  wave  height,  ,  over  the  roadway  (i.e., 

=  O.78  d^),  and  then  determining  the  fetch  required  to  generate  the 
height  of  the  breaking  wave.  The  fetch  necessary  to  generate  the  break¬ 
ing  wave  may  be  determined  from  the  appropriate  wave  forecasting  curves 
presented  in  the  SPM  (1977).  This  fetch  added  to  the  fetch  between  the 
roadway  and  the  project  site  can  then  be  used  to  estimate  the  wave 
height  and  wave  period  at  the  project  site. 

The  following  problem  will  illustrate  the  determinations  required 
to  estimate  the  wave  characteristics  at  a  project  site  for  the  case  when 
a  roadway  traversing  a  water  body  alters  the  effective  wave  fetch. 

EXAMPLE  PROBLEM 

GIVEN:  A  roadway  crosses  a  reservoir  midway  between  the  upwind  side  of 
the  reservoir  and  the  project  site  in  which  F  =  F  =  20,000  feet. 

The  fetch  F^  is  the  distance  between  the  roadway  and  the  upwind  side 
of  the  reservoir,  and  F^  is  the  downwind  fetch.  The  wind  speed,  U  , 
is  60  mph,  d^  =  d^  =  35  feet  and  the  still  water  level  is  5  feet  above 


20 


HTT80 


the  roadway.  At  the  project  site,  the  embankment  consists  of  a  rip- 
rapped  1  on  3  slope. 

FIND:  Determine  the  height  and  period  of  the  significant  wave  runup 
and  total  rise  in  water  level  at  the  project  site. 

SOLUTION: 

=  0.78(3^  =  0.78(5)  =  3-9  feet 

From  Figure  3-27,  Volume  I,  SPM  (1977) »  it  is  seen  that  the  fetch 
required  to  generate  the  3-9  foot  wave  is  about  12,000  feet;  therefore. 


F  =  20,000  +  12,000  =  32,000  feet 

For  the  modified  fetch  length,  Figure  3-27  reveals  that  Hg  =  5*5  feet 
and  T  =  U.7  seconds. 

The  deepwater  wave  length  according  to  equation  (5)  is 
Lo  =  5-12  T2  =  5-12  ( U . 7 ) 2  =  113  feet 
and  the  significant  wave  runup  (equation  17)  is 


R  = 


Hs 


5-5 


.l^Hs/Lo)1'2  cot  6  0.1  +(5-5/113)1/2  (3) 


=5.18  feet 


The  wind  setup  according  to  equation  (l8),  assuming  D  =  d  ,  is 


S  = 


* u  =  ( 60)2 ( 2 ) ( 32.000)  _ 


lloo  D  lloo  (35)  5280 


=0.89  feet 
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The  total  rise  in  water  level  at  the  project  site  is 
S  +  R  =  0.89  +  5.18  «  6.1  feet 

Finally,  it  should  be  noted  that  the  model  does  not  account  for 
wave  modification  as  a  result  of  refraction,  diffraction,  and  reflection. 
Reference  is  made  to  the  Shore  Protection  Manual  (1977)  for  treating 
such  wave  phenomena. 
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PART  II:  COMPUTER  FUNCTIONAL  DESCRIPTION 

1.  REVISION  LOG:  Program  coding  revised  June  78  to  accommodate  the 
CORPS  time-share  features. 


2.  FUNCTIONAL  FLOW  CHART:  i 


H7T0O 


HT7 


Pn 


Read:  Job  title,  site  ID,  number  radials,  type  of  slope 
cover,  use  overland  wind  speed  or  1  and  60  minute  wind  speed,) 
design  wave,  significant  wave  or  design  wave  is  other 
related  wave,  do  or  do  not  print  deep  water  wave  character¬ 
istics,  depth  at  toe  of  structure,  depth  along  fetch,  angle 
between  embankment  and  horizontal,  factor  for  reducing  wave 
runup  on  slightly  rougher  than  smooth  slope,  1  and  60  min 
wind  speed  or  overland  wind  speed  and  radials 


/  NSTART  X 

\L> 

S'  any  data  ^ 

\  T 

k1  r 

\  =  2  / 

V  before  run  , 

LQZ  =  2 

2h 
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H77B0I 


ENTRY  H77BO 


=  0  and  LQZ 


desired 


F  / 


READ:  all  inputs  -  job  title,  site  ID,  number  of 
radials  and  radials,  use  design  overland  wind 
speed  or  1  and  60  minute  wind  speed,  design  wave, 
significant  wave  or  design  wave  is  other  related 
wave,  do  or  do  not  print  deep  water. wave  char¬ 
acteristics,  depth  at  toe  structure,  depth  along 
fetch,  angle  between  embankment  and  horizontal, 
factor  for  reducing  wave  runup  on  slightly  rougher 
than  smooth  slope,  type  of  slope  cover,  and  save 
tabulated  output 


set  control 
SS7AKT  =  0 


/  save  ^ 
tabulated 
S'XPUtput/' 


set  control 
KFILE  =  2 

[  READ:  file  name 

NSAVE  =  1 

i  i 

IEZ13ZB3II 


set  control 
KFILE  =  1 


return 

-v,  =  1/ 


H7780W 


'return  X  F 
1 


N5TAKT=  0 


WRITE  to  tvFILE:  program  name  H7T80,  .job  title, 
site  ID,  number  of  radials,  type  of  slope  cover, 
whether  to  use  overland  or  1  and  60  minute  wind 
speed,  design  wave  factor,  print  option  for  deep 
water  wave  conditions,  depth  at  too,  depth  along 
fetch,  angle  between  embankment  and  horizontal, 
slope  factor,  overland  wind  speed  or  1  and  60 
minute  wind  speed  and  radials 
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3.  EQUIPMENT  AMD  OPERATING  SYSTEM:  The  basic  program  was  developed  on 
CDC  7600,  Southwestern  Division,  Dallas,  TX.  The  revised  computer 
version  in  this  abstract  was  developed  on  a  G635  time-share  system  in 
which  input/output  equipment  consisted  of  a  Model  33  remote  teletype. 

It  is  now  available  on  the  WES  G635,  Vicksburg,  MS;  HIS  66/80,  Macon, 

GA;  and  Boeing  CDC,  Seattle,  WA. 
h .  INPUT  REQUIREMENTS :  The  required  inputs  are: 

a.  Read,  in  the  computational  subroutine  H7?8o,  from  a  user 
input  data  file. 

b.  Entered,  in  subroutine  H7780I,  via  the  user's  time-share 
terminal  device  in  free  field  format.  The  inputs  are  passed 
to  the  computational  subroutine  H77&0  via  the  COMMON  state¬ 
ment  . 

5.  SECONDARY  STORAGE  INPUT  FORMAT:  The  formats  for  the  user  input 
data  file,  whether  it  already  exists  or  is  created  during  a  run  of  H7780, 
are : 


a.  READ: 

10260  FORMAT  (UX,A8) 

10350  FORMAT  (UX,15Al4AX,15AU/UX,5I3AX,>*F7.2,lX,2F7.?,7(AX, 
8F7.2)) 

Refer  to  line  numbers  10250  and  10330-103^0  of  the  source  listing, 
page  6l. 


b .  WRITE : 

16180  FORMAT  (I3,"H7780",2(/I3,1X,15AM/I3,1X,5I3/I3,1X,14F7.2/I3, 
IX , 2F7 . 2 , 7 ( / 1 3 , IX , 8F7 -  2 ) ) 

Refer  to  line  numbers  16150-16170  of  the  source  listing,  page  67. 
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6.  INPUT  DATA  DESCRIPTION:  The  following  names  are  used  for  the  input 
variables  in  program  HTJ80. 

COTTH  -  cotangent,  angle  between  horizontal  and  embankment 
slope.  Example  -  slope  of  1  on  2,  COTTH  =  2.0. 

D  -  average  water  depth  along  wind  fetch,  ft 

D1  -  water  depth  at  tow  of  structure,  ft 

FACT  -  factor  for  adjusting  slope  roughness  for  relatively 

smooth  slopes.  FACT  =  1.0  for  smooth  slopes;  =  0.85 
to  0.90  for  turfed  slopes 

FILEK(l)  -  8  character  variable  for  name  of  user  input  data  file 

FILEK(2)  -  8  character  variable  for  name  of  user  output  data  file 
for  graphics  and/or  other  use 

IDES  -  factor  lor  design  wave;  0,  design  wave  is  significant 

wave;  1,  design  wave  is  some  other  wave  in  a  representa¬ 
tive  spectrum 

I  MM  -  number  of  radials  in  evaluation  of  critical  effective- 

fetch;  15  minimum 

X{l)  -  length  of  radials  for  1  =  1,...,IMM;  miles 

ITAB  -  print  option  for  deepwater  wave  characteristics  of 

various  wind  speeds.  ITAB  =  1,  print  table  of  values; 

=  0,  no  print 

ITYPE  -  type  of  embankment  slope  cover.  ITY1E  =  -1,  riprap; 

=  0,  smooth;  =  1,  other  types 

IWIND  -  factor  for  design  wind.  1WINU  =0,  input  overland 

wind  speed;  =1,  input  1  and  60  minute  wind  speeds  as 
provided  in  ETL  1110-2-221 

UL  -  overland  wind  speed,  mph 

U2  -  maximum  1  minute  wind  speed,  mph 

U6l  -  maximum  60  minute  wind  speed,  mph 
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SITE  -  15  dimension,  U  character  variable  for  site 

identification 

TITLE  -  15  dimension,  1  character  variable  for  job  title 

7.  OUTPUT  DATA  DESCRIPTION :  The  following  names  are  used  for  the  out¬ 

put  variables  in  program  H7780. 

DL  -  design  deepwater  wave  length,  ft 

FE  -  critical  effective  wave  fetch,  miles 

FS  -  fetch  for  wind  setup,  miles 

H  -  significant  wave  height  (deepwater  characteristics),  ft 

HD  -  design  wave  height,  ft 

HP  -  wave  height  for  waves  other  than  significant  wave,  ft 

PPEP  -  wave  exceedance.  Exceedance  refers  to  the  percent  of 

wave  in  a  wave  spectrum  that  exceeds  a  given  value, 
percent 

RP  -  wave  runup  for  wave  other  than  significant  wave,  ft 

RS  -  significant  wave  runup,  ft  above  SWL 

RSS  -  total  increase  in  water  level,  ft  above  SWL 

S  -  wind  setup,  ft 

T  -  significant  wave  period  (deepwater  characteristics), 

sec 

TD  -  design  wave  period,  sec 

TDUR  -  wind  duration  (deepwater  characteristics),  min 

TDURD  -  design  wind  duration,  min 

UD  -  design  wind  speed,  mph 

-  over  water  wind  speed  (deepwater  characteristics), 
mph 


UF 


8.  PROGRAM  ERROR  MESSAGES: 


Messages  from  subroutine  TACHFILE.  TACHFILE  is  a  file  handling 
routine  for  CORPS  H  files. 

a.  If  file  name  INDATA  is  entered  as  an  existing,  but  does  not 
exist,  then 


NO  SUCH  FILE  INDATA  RE-ENTER  NAME,  Y  OR  N 

is  printed.  If  answer  is  Y,  then  the  new  file  name  is  read  and  the  file 
is  attached.  A  reply  of  N  returns  control  to  the  main  program. 

b.  If  INDATA  does  exist,  but  cannot  be  attached,  then 

FILE  PROB.  CALL  M  T  HEBLER  AT  88-5U2-2U03  AND  GIVE  THIS 
NUMBER  ISTAT  THANKS 

where  ISTAT  is  a  12  digit  octal  number.  Control  is  returned  to  the 
main  program. 

c.  If  INDATA  is  busy,  then 

FILE  INDATA  BUSY.  ANOTHER  FILE,  Y  OR  N 

is  printed.  If  answer  is  Y,  then  the  new  file  name  is  read  and  the  file 
is  attached.  If  N,  control  is  returned  to  the  main  program. 

d.  If  IN :DATA,  which  has  the  illegal  character  :  in  its  file  name 
occurs,  then 


ILLEGAL  CHAR  IN  FILE  NAME  IN: DATA 


RE-ENTER  Y  OR  N 
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is  printed.  If  answer  is  Y,  a  new  file  name  is  read  and  the  file 
attached.  If  N,  control  is  returned  to  the  main  program. 

e.  If  file  INDATA  is  an  existing  attached  file,  but  is  not  an 
input  file  for  H7780,  then 

FILE  INDATA  IS  NOT  AN  INPUT  FILE  FOR  H7780.  RE-ENTER 
Y  OR  N 

is  printed.  If  answer  is  Y,  then  the  new  file  is  attached  and  tested. 

If  N,  control  is  returned  to  the  main  program. 

f.  If  file  INDATA  is  an  existing  attached  file,  but  is  not  an 
output  file  for  H7780,  then 

FILE  INDATA  IS  NOT  AN  OUTPUT  FILE  FOR  H7780.  RE-ENTER 
Y  OR  N 

is  printed.  If  answer  is  Y,  the  new  file  is  attached  and  tested.  If  N, 
control  is  returned  to  the  main  program. 

9.  VARIABLE  DEFINITIONS: 
a.  Main  Program: 

HFILE  -  five  character  name  of  program;  passed  to  WESI.IB  count 
routine  HACCT 

LQZ  -  -  equal  1,  call  subroutine  H7780I  and  execute  all  input 

cues  and  reads;  =  2,  call  WESLIB  routine  RERUN  and 
enter  only  desired  inputs 

LQX  -  equal  1,  print  instructions  for  RERUN;  =  3,  no  print 

ZZZZZ  -  2  character;  =  RE,  rerun;  =  ST,  stop 
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b.  Subroutine  H7780: 


ADIFF  -  absolute  value  of  the  difference  between  the  trial 
transitional  depth  wave  length  and  the  value  of 
equation  (9),  page  13,  ft 

AI  -  wave  exceedance  in  percent 


ALEN 


ALI'J 


-  trial  transitional  depth  wave  length,  ft 

-  working  storage;  equal  to  the  expression 
equation  ( 4 )  page  11. 


of 


ALPHA ( I ) 


-  angle  between  a  specific  radial  line  and  the  central 
radial,  for  I  =  1,...,15,  radians 


AU  -  addition  counter;  adds  1  to  over  water  wind  speed  for 

for  each  print,  until  maximum  number  of  prints  reached 
for  wind  data  and  deepwater  wave  characteristics 

AU360  -  natural  log  of  wind  speed  for  duration  of  360  minutes, 

mph 

AU60  -  natural  log  of  maximum  60  minute  wind  speed,  mph 


AO 


Al 


a6o 

A360 

DLEN 

BLN 

C 


-  working  storage;  equal  to  the  expression 

exp(ln  Id  -  A  In  60)  of  equation  (13-4),  page  B-3, 
Appendix  B 


-  absolute  value  of  the  reciprocal  of  the  difference 
ratio  between  the  natural  logs  of  60  and  360  minute 
wind  speeds  to  the  natural  logs  of  60  and  360 


-  natural  log  of  6 0 

-  natural  log  of  360 


-  value  of  equation  (9)  page  13  at  trial  transitional 
depth  wave  length,  ft 

-  working  storage;  equal  to  the  expression  B  lnd^-j  of 

equation  (4)  page  11  \U  / 


-  equal  to  22/19;  changes  mph  to  fps 


COTTH  -  cotangent;  angle  between  embankment  slope  and  hori¬ 
zontal.  Example,  slope  of  1  on  2;  COTTH  =  2 
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D 

DIFF 

DL 

DLN 

DS 

D1 

F 

FACT 

FE 
FES 
FI  ( I ) 

FILEK(l) 

FILEK(2) 

FILET 

FS 

GD(I) 


-  average  water  depth  along  wind  fetch,  ft 

-  difference  between  trial  transitional  depth  wave 
length  and  the  value  of  equation  (9)  page  13  at  trial 
transitional  wave  length,  ft 

-  design  deepwater  wave  length,  ft 

-  working  storage;  equal  to  the  expression  D 
equation  (U)  page  11. 

-  depth  at  toe  of  structure  plus  wind  setup, 

-  depth  of  toe  of  structure,  ft 

-  effective  fetch,  miles 

-  factor  for  adjusting  slope  roughness  for  relatively 
smooth  slopes.  FACT  =  1,  smooth  slope;  =  .85  to  .9, 
turfed  slopes 

-  effective  fetch  (critical),  miles 

-  effective  fetch  (critical),  ft 

-  fetch  increments  used  in  determining  the  factor  for 
converting  overland  wind  speed  to  overwater  wind  speed 
for  I  =  1 , . . . ,7 

I  FI  I  FI  I  FI 

1  0  2  .5  3  1 

1*  2  5  3  6  It 

7  5 

-  8  character  name  of  input  data  file 

-  8  character  name  of  output  file  for  graphics  and/or 
other  use 

-  5  character  name  of  CORPS  H-file  used  to  build  input 
data  file  or  output  file  for  graphics  and/or  other  use 

-  fetch  for  wind  setup,  miles 

-  constants  for  equation  (M  page  11  for  I  =  1,...,5 
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GDU 

GFU 


GH 


GRAV 

GT 


GTA 

GTB 

G1 


G2 


G3 

GU 


I 

1 

u 


GD 

6.5882 

2.202U 


1 

2 
5 


GD 

0.0161 

0.8798 


GD 

0.3692 


fid 


working  storage;  equal  to  the  expression  of  equa¬ 
tions  (6)  and  (7)  pages  12  and  13  ^ 


-  working  storage;  equal  to  the  expression  of  equa- 

U2 

tions  (2),  (3),  and  (b)  page  11  and  equations  (6)  and 
(7)  pages  12  and  13 


-  working  storage;  equal  to  the  expression  tanh 

.  U2\ 

of  equation  (2)  page  11  and  equation 


.0125 


ifi 

(6)  page  12 


acceleration  of  gravity,  32.2  ft/sec^ 


-  working  storage;  equal  to  the  expression  tanh 

,-25\ 

•  077(^1  Jof  equation  (3)  page  1.1  and  equation  (7) 


page  13 

-  working  storage;  equal  to  the  expression 
tion  (9)  page  13 


fitf. 


of  equa- 


-  working  storage;  equal  to  the  expression  2nd  of  equa¬ 
tion  (9)  page  1; 


-  working  storage;  equal  to  the  expression  tanh 


,  s  - 


C) 


of  equation  (6)  page  12 


-  working  storage;  equal  to  the  expression  tanh 

•  375  \ 


■833(Sf)‘ 


of  equation  (7)  page  13 


-  working  storage;  equal  to  the  right  hand-side  of 
equation  (6)  page  12 

-  working  storage;  equal  to  the  right  hand-side  of 
equation  (7)  page  13 


37 


H7780 


H  ( I ) 

HD 

HFILE 
HP  ( I ) 

IDES 


[FIN 


II 


IJK 


IK 


IL 


I  MM 


INC 


lAB 


significant  wave  height  for  wind  data  and  deepwater 
wave  characteristics,  I  =  ft 

-  design  wave  height,  ft 

-  5  character  name  of  program's  CORPS  H-file  name  (H7780) 

-  wave  heights  for  waves  other  than  the  significant  wave, 
for  I  =  1, . . . ,20; ft 

-  factor  for  design  wave;  IDES  =  0,  design  wave  is 
oo-gnif leant  wave;  =  1,  design  wave  is  some  other 
related  wave 

-  switch;  =  0,  read  angles  between  specific  radial 
lines  and  central  radial  line  from  1+2°  down  to  0°, 

=  1,  read  back  up  from  6°  to  h2° 

-  counter  in  determination  of  effective  fetch;  =  1  if 
15  radial  lengths  input;  if  >  15  radial  lengths 
entered,  increased  one  each  time  effective  fetch 
scheme  is  entered  until  =  number  of  radial  lengths 
input  -15+1 

-  counter  in  determination  of  effective  fetch;  =  15  if 
15  radial  length  input;  increased  one  each  time 
effective  fetch  scheme  is  entered  until  =  number  of 
radial  lengths  input 


number  of  first  radial  length  corresponding  to  criti¬ 
cal  effective  fetch  to  be  printed 

-  number  of  last  radial  length  corresponding  to  critical 
effective  fetch  to  be  printed;  =  IK  +  it 

-  number  of  radial  lengths  input 

-  counter  used  in  determination  of  deepwater  wave  charac¬ 
teristics;  increased  by  1  until  a  fetch  increment  used 
in  determining  the  factor  for  converting  overland  wind 
speed  to  overwater  wind  speed  is  found  >_  to  the  criti¬ 
cal  effective  fetch 

-  option  for  printing  tabulation  of  deepwater  wave 
characteristics;  ITAB  =  0,  no  print;  =  1,  print 
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ITYPE  -  type  of  embankment  cover;  =  -1,  riprap;  =  0,  smooth; 

=  1,  other  types 

IUD  -  integer  used  to  convert  design  wind  speed  to  an 

integer  value,  mph 

IW1ND  -  factor  for  design  wind;  =  0,  input  overland  design 
wind;  =  1,  wind  values  input  from  data  provided  in 
ETL  1110-2-221  (l  and  60  minute  wind  speeds) 

IX  -  total  number  of  values  of  outputs  to  be  r?rinted  for 

wind  data  and  deepwater  wave  characteristics 

J  -  counter  in  determination  of  effective  fetch;  keeps 

track  of  which  angle  between  a  specific  radial  line 
and  the  central  radial  line  to  use  in  the  computations 

KFILE  -  logical  file  designator  for  input  data  file  and  output 
file  for  graphics  and/or  other  use;  =  1,  input  file; 

=  2,  output  file 

LQX  -  equal  1,  print  instruction  for  RERUN;  =  3,  no  print. 

No  function  in  this  subroutine  (H7?8o),  just  uses  COMMON 
to  be  passed  to  input  subroutine  H7780I 

LQZ  -  set  equal  to  1  if  no  permanent  input  data  file  exists; 

set  equal  to  2  if  file  exists.  If  =  1,  execute  cue 

"CHANGE  ANY  DATA  BEFORE  RUN,  Y  OR  N" 

NSTART  -  set  equal  to  1  if  permanent  data  file  exists  and  this 
is  first  entry  to  file;  set  equal  0  if  no  permanent 
file  or  permanent  file  has  already  been  entered. 

If  =  1,  execute  cue  "CHANGE  ANY  DATA  BEFORE  RUN,  Y 
OR  N" 

NSAVE  -  equal  1  or  2,  save  tabulated  output  for  graphics 
and/or  other  use;  =  0,  no  save 

P  -  wave  exceedance  in  hundredths 

PI  -  constant  =  3-1^159265 

PPEP(l)  -  wave  exceedance  in  percent  for  I  =  1,...,20 

R  -  ratio  of  overwater  wind  speed  to  overland  wind  speed 
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RAD 

RG 

RMAX 
RP(  I ) 

RS 

RSS 

S 

SCOSAL 

SINB 

SITE(I) 

SXA 

T(  I ) 

TD 

TDC 

TDF 

TDUR(I) 

TDURD 
TITLE (I) 


-  constant  =  tt/180;  converts  degrees  to  radians 

-  proportional  distance  of  critical  effective  fetch 
between  two  fetch  increments  Fl(l).  Example,  if 
critical  effective  fetch  is  3-75  and  Fl(l)  =  3  and 
FI(2)  =  4,  then  RG  =  .75 

-  greatest  radial  length,  miles 

-  wave  runup  for  waves  other  .nan  the  significant  wave, 
for  I  =  1, . . . ,20;  ft 

-  significant  wave  runup,  ft  above  SWL 

-  total  increase  in  water  level,  ft  above  SWL 

-  wind  setup,  ft 

-  sum  of  the  cosine  of  the  angles  between  a  specific 
radial  line  and  the  central  radial  line 

-  sine;  angle  between  embankment  slope  and  horizontal 

-  4  character  variable,  site  ID  name,  for  I  =  1,...,15. 
Can  have  total  of  60  characters  for  name 

-  working  storage;  used  in  determination  of  effective 
fetch,  miles 

-  significant  wave  period  for  deepwater  wave  character¬ 
istics,  for  I  =  1,...,IX,  sec 

-  design  wave  period,  sec 

-  wind  duration  as  expressed  by  equations  (B-2),  (B-3), 
and  (B-4),  pages  B-2  and  B-3,  Appendix  B 

-  value  of  equation  (4)  page  11 

-  wind  duration,  wind  data,  and  deepwater  wave 
characteristics  I  =  1,...,IX;  min 

-  design  wind  duration,  min 

-  4  character  variable,  job  title,  for  1  =  1,...,15. 

Can  have  total  of  60  characters  for  title 


Uo 


'-V, 
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U 

UD 

UF(l ) 

UL 

US 

U1 

U2 

U30 

U360 

U60 

u6i 

WR(l) 


X(I) 

XA 

ZA 


ZB 


-  wind  speed,  mph 

-  design  wind  speed,  mph 

-  wind  speed  overwater  for  wind  data  and  deepwater 
wave  characteristics,  1=1,..., IX;  mph 

-  input  overland  wind  speed,  mph 

-  wind  speed,  fps 


-  maximum  1  minute  wind  speed,  mph 

-  input  maximum  1  minute  wind  speed,  mph 

-  30  minute  wind  speed,  mph 

-  360  minute  wind  speed,  mph 

-  maximum  60  minute  wind  speed,  mph 

-  input  maximum  60  minute  wind  speed,  mph 

-  wind  factors  corresponding  to  fetch  increments  FI (I 


I  = 

I 

WR 

I 

WR 

I 

WR 

1 

1 

2 

1.08 

3 

1.13 

4 

7 

1.21 

1.30 

5 

1.26 

6 

1.28 

-  radial  lengths,  miles 

-  working  storage  in  determination  of  effective  fetch 
miles 


-  working  storage  used  in  determination  of  wave  runup: 

/H  \1/2 

equal  to  the  expression/  —  J  (1.58  -  2.35  tan  0) 
equation  (11)  page  14  \  s  / 


-  working  storage  used  in  determination  of  wave  runup 
equal  to  the  expression  .092  cot  0  -  .26  of  equa¬ 
tion  (11)  page  l4 


4l 
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ZC  -  working  storage  used  in  determination  of  wave  runup; 

equal  to  the  expression  sin  0  (5.95  tan  6  +1.5)  of 
equation  (10)  page  14 

ZE  -  value  of  equation  (ll)  page  14 

c.  Subroutine  H7780I:  Variables  are  as  explained  in  subroutine 
H7780,  except 

LQZ  -  equal  1,  execute  input  cues  and  reads;  =  2,  call  RERUN 

and  enter  only  desired  inputs 

JKL  -  direct  return  from  RERUN  to  desired  input  read 

KKK  -  total  number  of  inputs  passed  to  RERUN 

d.  Subroutine  H7780W:  Variables  are  as  explained  in  subroutine 
H7780  except 

NLIN  -  number  of  lines  in  data  file;  maximum  number  of  13 

LINE(l)  -  line  number  in  data  file;  start  at  100  and  incremented 
by  2  for  I  =  2, . . . ,NLIN 

M8  -  number  of  8  variable  lines  required  for  writing  radial 

length 

MR  -  number  of  radial  lengths  left  after  M8  lines  written 

e.  Subroutine  H7780H:  This  subroutine  writes  the  header  informa¬ 
tion  and  tabulated  output  to  the  output  file  for  graphics  and/or  other 
use.  The  header  information  may  be  seen  in  formats  20070  and  20100, 
line  numbers  20070-20080  and  20100-20140  of  the  source  listing,  page  68. 
Variables  are  as  explained  in  subroutine  H7780  except 

AFILE  -  10  character  variable  used  in  attaching  existing  input 
data  file 

LINE(l)  -  line  number  in  data  file;  start  at  100  and  increment 
by  2  for  I  =  1,...,IX;  maximum  number  50 
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10.  EXAMPLE  CASE:  The  example  problems  presented  were  taken  from  those 
given  in  ETL  1110-2-221  and  ETL  1110-2-8  (references  1  and  10 )  for 
Dension  and  McGee  Bend  Reservoirs.  Fetch,  average  depth,  and  slope 
inclination  data  used  correspond  to  the  values  given  in  those  ref¬ 
erences.  Wind  data  used  are  based  on  the  overland  maximum  one-minute 
and  sixty-minute  wind  speeds  obtained  from  the  appropriate  figures  in 
reference  1.  For  demonstration  purposes  only,  the  depths  at  the  toe  of 
the  structures  and  type  of  embankment  slopes  were  arbitrarily  selected. 
Examples  1  and  2  are  for  smooth  embankment  slopes  and  example  3  is  a  run 
with  different  data  files  using  the  RERUN  option  for  turfed  slopes, 
small  stepped  slopes,  and  riprapped  slopes. 

Example  1  Runs  program  with  no  permanent  data  file 

PT’N  VE«LIP/COFPC/H77tf0,p 
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INPUT  H7780  -  WAVE  RUNUP  AND  WIND  SETUP, COMPUTATIONAL  MODEL 

PERM  DATA  FILE; Y  OR  N 
=N 

WE  WILL  HELP  YOU  SETUP  YOUR  DATAFILE.  ENTER  FILE  NAME 
=  DH7  7  8  0 

AA-ENTER  JOB  TITLE  <  OR  =  60  CHARACTERS 
-DENSION  RESERVOIR 

AB-ENTER  SITE  ID  <  OR  =  60  CHARACTERS 
=STATION  A 

AC-ENTER  THE  NUMBER  OF  RADI ALS  USED  TO  DETERMINE  EFFECTIVE 
FETCH  DISTANCE; A  MINIMUM  OF  15  RADIALS  IS  REQUIRED 
=  15 

AD-ENTER  THE  15  RADIAL  LENGTHS  IN  MI L ES , SP EAR ATE  BY  COMMAS 

=2.23,2.21,2.92,3.2,5.05,4.58,5.45,8.02,7.86,7.53,2.11,1.71,1.24 

=1.25,1.22 

AE-USE  DESIGN  OVERLAND  WIND  SPEED  OR  1  MIN  AND  60  MIN  WIND  SPEED 
AS  IN  ETL  1110-2-221;0, DESIGN; 1 , 1  MIN  AND  60  MIN 
=  1 

AF-ENTER  1  MIN  AND  60  MIN  WIND  SPEEDS, MPH 
=65,45 

AG-ENTER  FACTOR  FOR  DESIGN  WAVE; 0 , DESIGN  WAVE  IS  SIGNIFICANT 
WAVE;1, DESIGN  WAVE  IS  OTHER  RELATED  WAVE 
=  1 

AH-ARE  DEEP  WATER  WAVE  CHARACTERISTICS  TO  BE  PRINTED;0,NO  PRINT ; 1 , PRINT 
=  1 

AI-ENTER  WATER  DEPTH  AT  TOE  OF  STRUCTURE, FT 
=  20 

AJ-ENTER  AVERAGE  WATER  DEPTH  ALONG  WIND  FETCH, FT 
=  50 

AK-ENTER  COTANGENT  OF  ANGLE  BETWEEN  EMBANKMENT  SLOPE  AND  HORIZONTAL 
=  2.5 

AL-ENTER  FACTOR  FOR  REDUCING  WAVE  RUNUP  ON  SMOOTH  SLOPE  WHERE  SLOPE  11 
SLIGHTLY  ROUGHER  THAN  SMOOTH 
=  1 

AM-ENTER  SLOPE  COVER  ;- 1 , RIPRAP ; 0 , SMOOTH ; 1 , OTHER 
=  0 

AN-SAVE  TABULATED  OUTPUT  FOR  GRAPHICS  OR  OTHER  USE 
=  Y 

FILE  NAME 
=PH7780 

CHANGE  ANY  DATA  BEFORE  RUN.Y  OR  N 
=  N 


OUTPUT  H7780  -  WAVE  RUNUP  AND  WIND  SETUP , COMPUT AT  I ONAL  MODEL 

DENSION  RESERVOIR 
STATION  A 


BASIC  DATA 

DEPTH  AT  TOE  OF  EMBANKMENT  =  20.00  FT 

AVERAGE  WATER  DEPTH  =  50.00  FT 

EMBANKMENT  SLOPE  =  2.50H:1V 

SMOuTH  SLOPE  MODIFICATION  FACTOR  =  1.00 

ONE  MINUTE  WIND  SPEED  =  65.00  MPH 

SIXTY  MINUTE  WIND  SPEED  =  45.00  MPH 


RADIALS  CORRESPONDING  TO  CRITICAL  EFFECTIVE  FETCH  IN  MILES 

2.23  2.21  2.92  3.20  5.05  4.58  5.45  8.02 

7.86  7.53  2.11  1.71  1.24  1.25  1.22 
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COMPUTATIONAL  RESULTS 
EFFECTIVE  WAVE  FETCH  =  3.75  MILES 


TRANSITIONAL  WATER  WAVE  CONDITIONS  USED  FOR  PRESENT  PROBLEM 

EMBANKMENT  SLOPE  -  SMOOTH 

DESIGN  WIND  SPEED  =  63.00  MPH 

DESIGN  WIND  DURATION  =  36.34  MIN 

DESIGN  WAVE  HEIGHT  =  5.29  FT 

DESIGN  WAVE  PERIOD  =  4.47  SEC 

DESIGN  DEEP  WATER  WAVE  LENGTH  =  102.80  FT 

FETCH  FOR  WIND  SETUP  =  7.50  MILES 

WIND  SETUP  =  0.42  FT 

SIGNIFICANT  WAVE  RUNUP  =  9.86  FT  ABOVE  SWL 

TOTAL  INCREASE  IN  WATER  LEVEL  =  10.29  FT  ABOVE  SWL 


WIND  DATA  AND  DEEP  WATER  WAVE  CHARACTERISTICS 

OVER  WATER  WIND  SIGNIFICANT  SIGNIFICANT 

WIND  SPEED  DURATION  WAVE  HEIGHT  WAVE  PERIOD 

(MPH)  (MIN)  (FT)  (SEC) 

50.00 
51.00 
52.00 
53.00 
54.00 
55.00 
56 . 00 
57.00 
58.00 
59.00 
60.00 
61.00 
62.00 
63.00 
64.00 
65.00 
66.00 
67.00 
68.00 
69.00 
70.00 
71.00 
72.00 
73.00 
74 .00 
75.00 
76.00 
77.00 
78.00 
79.00 
80.00 
81.00 


>*5 


40.97 

4.38 

4.27 

40.55 

4.48 

4.32 

40.15 

4.58 

4.36 

39.75 

4.68 

4.40 

39.37 

4.79 

4.45 

39.00 

4.89 

4.49 

38.63 

4.99 

4.53 

38 . 28 

5.10 

4.57 

37 . 94 

5.20 

4.61 

37.60 

5.31 

4.65 

37.27 

5.41 

4.69 

36 . 96 

5.52 

4.74 

36.65 

5.62 

4.78 

36 .34 

5.73 

4.81 

36.05 

5.83 

4.85 

35.76 

5.94 

4.89 

35.48 

6  .  05 

4.93 

35.20 

6.15 

4 . 97 

34.93 

6.26 

5.01 

34.67 

6  .  37 

5.05 

34.41 

6 . 47 

5.08 

34.16 

6.58 

5.12 

33.92 

6.69 

5.16 

33.68 

6.80 

5.19 

33.44 

6  .91 

5.23 

33.21 

7.01 

5.27 

32.93 

7.12 

5.30 

32.76 

7.23 

5.34 

32.54 

7  .  34 

5.37 

32.33 

7.45 

5.41 

32.12 

7.56 

5.44 

31 . 91 

7.67 

5.48 
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WAVE  HEIGHTS  AND  WAVE  RUNUP  FOR  WAVES  OTHER  THAN  THE  SIGNIFICANT  WAVE 
(NOTE-  THE  TERM  WAVE  EXCEEDENCE  REFERS  TO  THE  PERCENT  OF  WAVES  IN 
A  WAVE  SPECTRUM  THAT  EXCEEDS  A  GIVEN  VALUE.) 


WAVE 

WAVE 

WAVE 

EXCEEDANCE 

HEIGHT 

RUNUP 

IN  PERCENT 

(FT) 

(FT) 

1 .  00 

8 .03 

14.97 

2.00 

7.40 

13.80 

3.00 

7.01 

13.06 

4.00 

6.71 

12.51 

5.00 

6  .  48 

12.07 

6.00 

6.28 

11.70 

7.00 

6.10 

11.38 

8 . 00 

5.95 

11.09 

9.00 

5.81 

10.82 

10.00 

5.68 

10.58 

11.00 

5.56 

10.36 

12.00 

5.45 

10.16 

13.00 

5.34 

9  .  96 

14.00 

5.25 

9.78 

15.00 

5.15 

9.61 

16.00 

5.06 

9.44 

17.00 

4.98 

9.29 

18.00 

4 .90 

9.13 

19.00 

4.82 

8.99 

20.00 

4.75 

8.85 

ENTER  RERUN  OR  STOP 
=  STOP 


The  user  named  input  data  file  DH7780  is  now  saved  as  a  permanent  data 


on  the  user’s  ID.  The  following  is  a  list  of  DH7780. 


100 

H7780 

102 

DENS  ION 

RESERVOIR 

104 

STATION 

A 

106 

15  0 

1  1  1 

108 

20.00 

50.00  2.50 

1.00 

110 

65.00 

45.00 

112 

2.23 

2.21  2.92 

3.20 

114 

7.86 

7.53  2.11 

1.71 

5.05  4.58  5.45  8.02 

1.24  1.25  1.22 


The  file  PH7780,  a  user  named  output  file  for  graphics  and/or  other 


use,  is  also  saved  on  his/her  ID.  The  following  is  a  list  of  PH7780. 


be 
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10  H77SC  07  07  EDGE 

11  (07C/),<3X,7(1X.F7.2>>> 

12  32  32  32  32  20  20  20 

13  1  1  1  1  5  5  5 

14  CURVE  DESIGNATIONS  FOR  H7780  ARE-' 

15  1 =0V  ER  WATER  WIND  SPEED  2-WIND  DURATION 

16  3-SIGNIFICANT  WAVE  HEIGHT  4 =  S IGN I F I C ANT  WAVE  PERIOD 

17  5-WAVE  EXCEEDANCE  6=WAVE  HEIGHT 

18  7  =WAV  E  RUNUP 

19  UNITS  FOR  ABOVE  VARIABLES  ARE'- 


20 

MPH  =  1  MIN 

1  =  2  SEC 

=  4  FT  =  3 

,6,7  1 

5ERCENT  =  5 

100 

50.00 

40.97 

4 . 38 

4.27 

1  .  CO 

8.03 

14.97 

102 

51.00 

40 .55 

4.48 

4  .  32 

2.00 

7.40 

13.80 

104 

52.00 

40.15 

4  .  58 

4 . 36 

3.00 

7.01 

13.06 

106 

53.00 

39.75 

4 .68 

4.40 

4.00 

6.71 

12.51 

108 

54.00  ' 

39.37 

4.79 

4.45 

5.00 

6  .  48 

12.07 

110 

55.00 

39.00 

4.89 

4.49 

6.00 

6 . 28 

11.70 

112 

56.00 

38.63 

4.99 

4 .53 

7.00 

6.10 

11.38 

114 

57.00 

33.28 

5.10 

4 . 57 

8.00 

5.95 

11.09 

116 

58.00 

3  7.94 

5 .20 

4.61 

9.00 

5.81 

10.82 

118 

59.00 

37.60 

5.31 

4.65 

10.00 

5 . 68 

10.58 

120 

60.00 

37.27 

5.41 

4.69 

11.00 

5.56 

10.36 

122 

61.00 

36  .  96 

5.52 

4.74 

12.00 

5.45 

10.16 

124 

62.00 

36.65 

5.62 

4 . 78 

13.00 

5.34 

9.  96 

126 

63.00 

36 . 34 

5.73 

4.81 

14.00 

5.25 

9.78 

128 

64.00 

36  .05 

5.83 

4.85 

15.00 

5.15 

9.61 

130 

65.00 

35.76 

5.94 

4.89 

16.00 

6  .  06 

9.44 

132 

66.00 

35.48 

6.05 

4  .  93 

17.00 

4  .  93 

9.29 

134 

67  .00 

35.20 

6.15 

4 . 97 

18.00 

4  .  90 

9.13 

1  36 

68.00 

34 . 93 

6.26 

5.01 

19.00 

4.82 

8 . 99 

138 

69.00 

34.67 

6  .  37 

5.05 

20.00 

4 .75 

8.85 

140 

70.00 

34.41 

6.47 

5  .  08 

0  . 

0. 

0  . 

142 

71.00 

34.16 

6  .  58 

5.12 

0  . 

0  . 

0  . 

144 

72.00 

33.9? 

6.69 

5.16 

0  . 

0  . 

0  . 

146 

73.00 

33.68 

6.80 

5.19 

0  . 

0  . 

0  . 

1  48 

74.00 

3  3.44 

6.91 

5.2  5 

0  . 

0  . 

0  . 

150 

75.00 

33.21 

7.01 

5.27 

0  . 

0  . 

0  . 

152 

76.00 

32 . 98 

7.12 

5  .  30 

0  . 

0  . 

0  . 

154 

77.00 

32.76 

7.23 

5.34 

0  . 

0  . 

0  . 

156 

78.00 

32.54 

7  .  34 

5.37 

0  . 

0  . 

0  . 

158 

79.00 

32.33 

7.45 

5.41 

0  . 

0  . 

0  . 

160 

80.00 

32.12 

7  .  56 

5 .44 

0  . 

0  . 

0  . 

162 

81.00 

31.91 

7.67 

5.48 

0  . 

0  . 

0  . 

Example  ?  Run  program  u;; j rig  an  existing  permanent.  data  file 


>>7 


H778C 


INPUT  H7780  -  WAVE  RUNUP  AND  WIND  SETUP , COMPUTATIONAL  MODEL 

PERM  DATA  FILE;Y  OR  N 
=  Y 

DATAFILE  NAME 
-DH77801 

CHANGE  ANY  DATA  BEFORE  RUN , Y  OR  N 
-N 


OUTPUT  H7780  -  WAVE  RUNUP  AND  WIND  SETUP , COMPUTATIONAL  MODEL 

MCGEE  BEND  RESERVOIR 
STATION  A 


BASIC  DATA 

DEPTH  AT  TOE  OF  EMBANKMENT  =  20.00  FT 

AVERAGE  WATER  DEPTH  =  73.00  FT 

EMBANKMENT  SLOPE  =  2.50HUV 

SMOOTH  SLOPE  MODIFICATION  FACTOR  =  1.00 

ONE  MINUTE  WIND  SPEED  =  55.00  MPH 

SIXTY  MINUTE  WIND  SPEED  =  50.00  MPH 


RADI ALS  CORRESPONDING  TO  CRITICAL  EFFECTIVE  FETCH  IN  MILES 

1.33  2.70  2.89  3.22  6.49  12.64  16.95  16.71 

12.31  11.74  6.39  5.87  5.59  3.46  3.41 


COMPUTATIONAL  RESULTS 
EFFECTIVE  WAVE  FETCH  =  7.49  MILES 


TRANSITIONAL  WATER  WAVE  CONDITIONS  USED  FOR  PRESENT  PROBLEM 

EMBANKMENT  SLOPE  -  SMOOTH 

DESIGN  WIND  SPEED  =  64.00  MPH 

DESIGN  WIND  DURATION  =  61.01  MIN 

DESIGN  WAVE  HEIGHT  =  7.18  FT 

DESIGN  WAVE  PERIOD  =  5.31  SEC 

DESIGN  DEEP  WATER  WAVE  LENGTH  =  144.43  FT 

FETCH  FOR  WIND  SETUP  =  14.98  MILES 

WIND  SETUP  -  0 . 60  FT 

SIGNIFICANT  WAVE  RUNUP  =  14.18  FT  ABOVE  SW1 

TOTAL  INCREASE  IN  WATER  LEVEL  =  14.78  FT  ABOVE  SWL 


UR 


t 


H7780 


WIND  DATA  AND  DEEP  WATER  WAVE  CHARACTERISTICS 


OVER  WATER 

WIND 

SIGNIFICANT 

SIGNIFICANT 

WIND  SPEED 

DURATION 

WAVE  HEIGHT 

WAVE  PERIOD 

CMPM  J 

(MIN) 

(FT) 

(SEC) 

57  .00 

64.80 

6.81 

5.40 

58.00 

64.22 

6  .  95 

5.44 

59.00 

63.65 

7.09 

5.49 

60.00 

63.09 

7.23 

5 . 54 

61.00 

62.55 

7.37 

5.59 

62.00 

62.03 

7 .51 

5.64 

63.00 

61.51 

7.65 

5.69 

64 .00 

61.01 

7.79 

5.73 

65.00 

60.52 

7  .  93 

5.78 

66.00 

60.04 

8.08 

5.83 

67.00 

59.57 

8.22 

5.87 

68.00 

59.11 

8 . 36 

5.92 

69.00 

58.67 

8.51 

5.96 

70.00 

58.23 

8.65 

6.01 

WAVE  HEIGHTS  AND  WAVE  RUNUP  FOR  WAVES  OTHER  THAN  THE  SIGNIFICANT  WAVE 
(NOTE-  THE  TERM  WAVE  EXCEEDENCE  REFERS  TO  THE  PERCENT  OF  WAVES  IN 
A  WAVE  SPECTRUM  THAT  EXCEEDS  A  GIVEN  VALUE.) 


WAVE 

WAVE 

WAVE 

EXCEEDANCE 

HEIGHT 

RUNUP 

IN  PERCENT 

(FT) 

(FT) 

1.00 

10.90 

21.52 

2.00 

10.05 

19.84 

3.00 

9.51 

18.78 

4.00 

9.11 

17.99 

5-00 

8.79 

17.36 

6.00 

8 . 52 

16.82 

7  .00 

8.28 

16  .  35 

8.00 

8.07 

15.94 

9.00 

7.88 

15.56 

10.00 

7.71 

15.22 

11.00 

7.55 

14.90 

12.00 

7.40 

14.60 

13.00 

7.25 

14.32 

14.00 

7.12 

14.06 

15.00 

7.00 

13.81 

16.00 

6.88 

13.58 

17.00 

6.76 

13.35 

18.00 

6.65 

13.13 

19.00 

6  .54 

12.92 

20 . 00 

6.44 

12.72 

ENTER  RERUN  OR  STOP 
-STOP 
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The  following  is  a  list  of  the  permanent  input  data  file  PK77801  used 
in  example  2. 


100  H7780 

102  MCGEE  BEND  RESERVOIR 
104  STATION  A 


106 

15  0 

1  1  1 

108 

20.00 

73.00 

2.50 

1.00 

no 

55.00 

50.00 

112 

1  .  33 

2.70 

2.89 

3 . 22 

6  .49 

12.64 

16.95  16 

114 

12.31 

11.74 

6  .39 

5.87 

5.59 

3.46 

3.41 

Example  3  Run  program  using  different  existing  data  files  and  the 


rerun  option. 


1 


INPUT  H7780  -  WAVE  RUNUP  AND  WIND  SETUP , COMPUTATIONAL  MODEL 


PERM  DATA  FILE; Y  OR  N 

=  Y 

DATAFILE  NAME 
-DH  7780 

CHANGE  ANY  DATA  BEFORE  RUN , Y  OR  N 
-  Y 

RERUN  OPTION  PERMITS  YOU  TO  CHANGE  ANY  OR  ALL  INPUT  VARIABLES. 
AT  >>>=  QUE  TYPE  IN  THE  TWO  L ETTERS ( AA , AB , ETC .  )  CORRESPONDING 
TO  THE  VARIABLES  YOU  WISH  TO  CHANGE.  THEN  AT  NEXT  -  QUO,  ENTER 
THE  NUMERICAL  VALUE.  TO  TERMINATE  DATA  ENTRY,  TYPE  A  CARRIAGE 
RETURN  AT  »>  =  QUE. 


J 


>;■> 


-AG 

AG-ENTER  FACTOR  FOR  DESIGN  WAVE; 0 , DESIGN  WAVE  IS  SIGNIFICANT 
WAVE;1, DESIGN  WAVE  IS  OTHER  RELATED  WAVE 
=  0 


>>> 


-AH 

AH-ARE  DEEP  WATER  WAVE  CHARACTERISTICS  TO  BE  PRINTED;!), NO  PRINT  ;  1  ,  PRINT 
=  0 
>>> 

-AL 

AL-ENTER  FACTOR  FOR  REDUCING  WAVE  RUNUP  ON  SMOOTH  SLOPE  WHERE  SLOPE  IS 
SLIGHTLY  ROUGHER  THAN  SMOOTH 
-  .  9 
>>> 


H773f 

OUTPUT  H7780  -  WAVE  RUNUP  AND  WIND  SETUP , COMPUTAT IONAL  MODEL 

PENSION  RESERVOIR 
STATION  A 


BASIC  DATA 

DEPTH  AT  TOE  OF  EMBANKMENT  =  20.00  FT 

AVERAGE  WATER  DEPTH  =  50.00  FT 

EMBANKMENT  SLOPE  -  2.50h':lV 

SMOOTH  SLOPE  MODIFICATION  FACTOR  =  0.90 

ONE  MINUTE  WIND  SPEED  =  65.00  MPH 

SIXTY  MINUTE  WIND  SPEED  =  45.00  MPH 


RADI ALS  CORRESPONDING  TO  CRITICAL  EFFECTIVE  FETCH  IN  MILES 


2.23  2.21 

7.36  7.53 


2.92  3.20 

2.11  1.71 


5.05  4.58 

1.24  1.25 


5.45  8.02 

1 .22 


COMPUTATIONAL  RESULTS 
EFFECTIVE  WAVE  FETCH  -  3.75  MILES 


TRANSITIONAL  WATER  WAVE  CONDITIONS  USED  FOR  PRESENT  PROBLEM 

EMBANKMENT  SLOPE  -  MODIFIED  SMOOTH 

DESIGN  WIND  SPEED  =  63.00  MPH 

DESIGN  WIND  DURATION  =  36.34  MIN 

DESIGN  WAVE  HEIGHT  -  5.29  FT 

DESIGN  WAVE  PERIOD  =  4.47  SEC 

DESIGN  DEEP  WATER  WAVE  LENGTH  =  102.80  FT 

FETCH  FOR  WIND  SETUP  =  7.50  MILES 

WIND  SETUP  -  0.42  FT 

SIGNIFICANT  WAVE  RUNUP  =  8.88  FT  ABOVE  SWL 

TOTAL  INCREASE  IN  WATER  LEVEL  =  9.30  FT  ABOVE  SWL 

ENTER  RERUN  OR  STOP 
-RERUN 

INPUT  H7780  -  WAVE  RUNUP  AND  WIND  SETUP , COMPUTATIONAL  MODEL 

CHANGE  SPECIFIC  ITEMS  ON  EXISTING  DATA  SET  OR  NEW  DATA  SET 
NEW  Y  OR  N 
=  Y 

PERM  DATA  FILE;Y  OR  N 
=  Y 

DATAFILE  NAME 
: DH  7  7  3  0 1 

CHANGE  ANY  DATA  BEFORE  RUN , Y  OR  N 
-Y 
>>> 

-AG 

AG-ENTER  FACTOR  FOR  DESIGN  WAV E ; 0 , DES I GN  WAVE  IS  SIGNIFICANT 
WAVE ; 1 , DESIGN  WAVE  IS  OTHER  RELATED  WAVE 
-0 
>>> 

-AH 

AH-ARE  DEEP  WATER  WAVE  CHARACTERISTICS  TO  BE  PRINTED;0,NO  PR  I N T ; 1 , PR  I  NT 
=  0 
>>> 

=  AL 
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AL-ENTER  FACTOR  FOR  REDUCING  WAVE  RUNUP  ON  SMOOTH  SLOPE  WHERE  SLOPE  IS 
SLIGHTLY  ROUGHER  THAN  SMOOTH 
=  .9 
>>> 


OUTPUT  H7780  -  WAVE  RUNUP  AND  WIND  SETUP , COMPUTATIONAL  MODEL 

MCGEE  BEND  RESERVOIR 
STATION  A 

BASIC  DATA 

DEPTH  AT  TOE  OF  EMBANKMENT  =  20.00  FT 

AVERAGE  WATER  DEPTH  =  73.00  FT 

EMBANKMENT  SLOPE  =  2.50HUV 

SMOOTH  SLOPE  MODIFICATION  FACTOR  =  0.90 

ONE  MINUTE  WIND  SPEED  =  55.00  MPH 

SIXTY  MINUTE  WIND  SPEED  =  50.00  MPH 


RADIALS  CORRESPONDING  TO  CRITICAL  EFFECTIVE  FETCH  IN  MILES 

1.33  2.70  2.89  3.22  6.49  12.64  16.95  16.71 

12.31  11.74  6.39  5.87  5.59  3.46  3.41 


COMPUTATIONAL  RESULTS 
EFFECTIVE  WAVE  FETCH  =  7.49  MILES 


TRANSITIONAL  WATER  WAVE  CONDITIONS  USED  FOR  PRESENT  PROBLEM 

EMBANKMENT  SLOPE  -  MODIFIED  SMOOTH 

DESIGN  WIND  SPEED  =  64.00  MPH 

DESIGN  WIND  DURATION  =  61.01  MIN 

DESIGN  WAVE  HEIGHT  =  7.18  FT 

DESIGN  WAVE  PERIOD  =  5.31  SEC 

DESIGN  DEEP  WATER  WAVE  LENGTH  =  144.43  FT 

FETCH  FOR  WIND  SETUP  =  14.98  MILES 

WIND  SETUP  =  0.60  FT 

SIGNIFICANT  WAVE  RUNUP  =  12.76  FT  ABOVE  SWL 

TOTAL  INCREASE  IN  WATER  LEVEL  =  13.36  FT  ABOVE  SWL 

ENTER  RERUN  OR  STOP 
=RERUN 

INPUT  H7780  -  WAVE  RUNUP  AND  WIND  S ETUP , COMPUTAT I ONAL  MODEL 

CHANGE  SPECIFIC  ITEMS  ON  EXISTING  DATA  SET  OR  NEW  DATA  SET 
NEW  Y  OR  N 
=  N 
>>> 

=  AL 

AL-ENTER  FACTOR  FOR  REDUCING  WAVE  RUNUP  ON  SMOOTH  SLOPE  WHERE  SLOPE  IS 
SLIGHTLY  ROUGHER  THAN  SMOOTH 
=  .75 
>>> 
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OUTPUT  H7780  -  WAVE  RUNUP  AND  WIND  SETUP , COMPUTATIONAL  MODEL 

MCGEE  BEND  RESERVOIR 
STATION  A 

BASIC  DATA 

DEPTH  AT  TOE  OF  EMBANKMENT  =  20.00  FT 

AVERAGE  WATER  DEPTH  =  73.00  FT 

EMBANKMENT  SLOPE  =  2.50H:1V 

SMOOTH  SLOPE  MODIFICATION  FACTOR  =  0.75 

QNE  MINUTE  WIND  SPEED  =  55.00  MPH 

SIXTY  MINUTE  WIND  SPEED  =  50.00  MPH 


RADI ALS  CORRESPONDING  TO  CRITICAL  EFFECTIVE  FETCH  IN  MILES 

1.33  2.70  2.89  3.22  6.49  12.64  16.95  16.71 

12.31  11.74  6.39  5.87  5.59  3.46  3.41 


COMPUTATIONAL  RESULTS 
EFFECTIVE  WAVE  FETCH  =  7.49  MILES 


TRANSITIONAL  WATER  WAVE  CONDITIONS  USED  FOR  PRESENT  PROBLEM 

EMBANKMENT  SLOPE  -  MODIFIED  SMOOTH 

DESIGN  WIND  SPEED  =  64.00  MPH 

DESIGN  WIND  DURATION  =  61.01  MIN 

DESIGN  WAVE  HEIGHT  =  7.18  FT 

DESIGN  WAVE  PERIOD  =  5.31  SEC 

DESIGN  DEEP  WATER  WAVE  LENGTH  =  144.43  FT 

FETCH  FOR  WIND  SETUP  =  14.98  MILES 

WIND  SETUP  -  0.60  FT 

SIGNIFICANT  WAVE  RUNUP  =  10.64  FT  ABOVE  SWL 

TOTAL  INCREASE  IN  WATER  LEVEL  =  11.24  FT  ABOVE  SWL 

ENTER  RERUN  OR  STOP 
=RERUN 

INPUT  H7780  -  WAVE  RUNUP  AND  WIND  S ETUP , COMPUTATIONAL  MODEL 

CHANGE  SPECIFIC  ITEMS  ON  EXISTING  DATA  SET  OR  NEW  DATA  SET 
NEW  Y  OR  N 
-N 
>>> 

=  AL 

AL-ENTER  FACTOR  FOR  REDUCING  WAVE  RUNUP  ON  SMOOTH  SLOPE  WHERE  SLOPE  IS 
SLIGHTLY  ROUGHER  THAN  SMOOTH 
=  1 
>>> 

=  AM 

AM-ENTER  SLOPE  COVER  ; -1 , RIPRAP ; 0 , SMOOTH ; 1 , OTHER 
=  -l 
>>> 
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OUTPUT  H7780  -  WAVE  RUNUP  AND  WIND  SETUP, COMPUTATIONAL  MODEL 

MCGEE  BEND  RESERVOIR 
STATION  A 


BASIC  DATA 

DEPTH  AT  TOE  OF  EMBANKMENT  =  20.00  FT 

AVERAGE  WATER  DEPTH  =  73.00  FT 

EMBANKMENT  SLOPE  =  2. 50H:1V 

SMOOTH  SLOPE  MODIFICATION  FACTOR  =  1.00 

ONE  MINUTE  WIND  SPEED  =  55.00  MPH 

SIXTY  MINUTE  WIND  SPEED  =  50.00  MPH 


RADIALS  CORRESPONDING  TO  CRITICAL  EFFECTIVE  FETCH  IN  MILES 

1.33  2.70  2.89  3.22  6.49  12.64  16.95  16.71 

12.31  11.74  6.39  5.87  5.59  3.46  3.41 


COMPUTATIONAL  RESULTS 
EFFECTIVE  WAVE  FETCH  =  7.49  MILES 


TRANSITIONAL  WATER  WAVE  CONDITIONS  USED  FOR  PRESENT  PROBLEM 
EMBANKMENT  SLOPE  -  RIPRAP 


DESIGN  WIND  SPEED  =  64.00  MPH 

DESIGN  WIND  DURATION  =  61.01  MIN 

DESIGN  WAVE  HEIGHT  =  7.18  FT 

DESIGN  WAVE  PERIOD  =  5.31  SEC 

DESIGN  DEEP  WATER  WAVE  LENGTH  =  144.43  FT 

FETCH  FOR  WIND  SETUP  =  14.98  MILES 

WIND  SETUP  =  0 .60  FT 

SIGNIFICANT  WAVE  RUNUP  =  7.50  FT  ABOVE  SWL 

TOTAL  INCREASE  IN  WATER  LEVEL  =  8.10  FT  ABOVE  SWL 

ENTER  RERUN  OR  STOP 
=RERUN 

INPUT  H7780  -  WAVE  RUNUP  AND  WIND  S ETUP , COMPUTATI ONAL  MODEL 

CHANGE  SPECIFIC  ITEMS  ON  EXISTING  DATA  SET  OR  NEW  DATA  SET 
NEW  Y  OR  N 
=Y 

PERM  DATA  FILE;Y  OR  N 
=  Y 

DATAFILE  NAME 
-DH7780 

CHANGE  ANY  DATA  BEFORE  RUN, Y  OR  N 
=Y 
>>> 

=  AL 

AL-ENTER  FACTOR  FOX  REDUCING  WAVE  RUNUP  ON  SMOOTH  SLOPE  WHERE  SLOPE  IS 
SLIGHTLY  ROUGHER  THAN  SMOOTH 
=  .75 
>>> 


HTT80 


OUTPUT  H7780  -  WAVE  RUNUP  AND  WIND  SETUP , COMPUTAT IONAL  MODEL 

DfcNS ION  RESERVOIR 
STATION  A 

BASIC  DATA 

DEPTH  AT  TOE  OF  EMBANKMENT  =  20.00  FT 

AVERAGE  WATER  DEPTH  =  50.00  FT 

EMBANKMENT  SLOPE  =  2.50H:1V 

SMOOTH  SLOPE  MODIFICATION  FACTOR  =  0.75 

ONE  MINUTE  WIND  SPEED  =  65.00  MPH 

SIXTY  MINUTE  WIND  SPEED  =  45.00  MPH 


RADIALS  CORRESPONDING  TO  CRITICAL  EFFECTIVE  FETCH  IN  MILES 

2.23  2.21  2.92  3.20  5.05  4.58  5.45  8.02 

7.86  7.53  2.11  1.71  1.24  1.25  1.22 


COMPUTATIONAL  RESULTS 
EFFECTIVE  WAVE  FETCH  =  3.75  MILES 


TRANSITIONAL  WATER  WAVE  CONDITIONS  USED  FOR  PRESENT  PROBLEM 

EMBANKMENT  SLOPE  -  MODIFIED  SMOOTH 

DESIGN  WIND  SPEED  =  63.00  MPH 

DESIGN  WIND  DURATION  =  36.34  MIN 

DESIGN  WAVE  HEIGHT  =  5.29  FT 

DESIGN  WAVE  PERIOD  =  4.47  SEC 

DESIGN  DEEP  WATER  WAVE  LENGTH  =  102.80  FT 

FETCH  FOR  WIND  SETUP  =  7.50  MILES 

WIND  SETUP  =  0 .42  FT 

SIGNIFICANT  WAVE  RUNUP  =  7.40  FT  ABOVE  SWL 

TOTAL  INCREASE  IN  WATER  LEVEL  =  7.82  FT  ABOVE  SWL 

ENTER  RERUN  OR  STOP 
=RERUN 

INPUT  H7780  -  WAVE  RUNUP  AND  WIND  SETUP , COMPUT AT  I  ON AL  MODEL 

CHANGE  SPECIFIC  ITEMS  ON  EXISTING  DATA  SET  OR  NEW  DATA  SET 
NEW  Y  OR  N 
=N 
>>> 

-AL 

AL-ENTER  FACTOR  FOR  REDUCING  WAVE  RUNUP  ON  SMOOTH  SLOPE  WHERE  SLOPE  IS 
SLIGHTLY  ROUGHER  THAN  SMOOTH 
=  1 
>>> 

=  AM 

AM-ENTER  SLOPE  COVER  ; -1 . RIPRAP ; 0 , SMOOTH ; 1 , OTHER 

=  -l 

>>> 
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OUTPUT  H7780  -  WAVE  RUNUP  AND  WIND  SETUP , COMPUTATIONAL  MODEL 

DENSION  RESERVOIR 
STATION  A 


BASIC  DATA 

DEPTH  AT  TOE  OF  EMBANKMENT  =  20.00  FT 

AVERAGE  WATER  DEPTH  =  50.00  FT 

EMBANKMENT  SLOPE  =  2.50H'-1V 

SMOOTH  SLOPE  MODIFICATION  FACTOR  =  1.00 

ONE  MINUTE  WIND  SPEED  =  65.00  MPH 

SIXTY  MINUTE  WIND  SPEED  =  45.00  MPH 


RADI ALS  CORRESPONDING  TO  CRITICAL  EFFECTIVE  FETCH  IN  MILES 


2.23 

2.21 

2.92 

3.20 

5.05 

4.58 

5.45 

8.02 

7.86 

7.53 

2.11 

1.71 

1.24 

1.25 

1.22 

COMPUTATIONAL  RESULTS 
EFFECTIVE  WAVE  FETCH  =  3.75  MILES 


TRANSITIONAL  WATER  WAVE  CONDITIONS  USED  FOR  PRESENT  PROBLEM 
EMBANKMENT  SLOPE  -  RIPRAP 


DESIGN  WIND  SPEED  =  63.00  MPH 

DESIGN  WIND  DURATION  =  36.34  MIN 

DESIGN  WAVE  HEIGHT  =  5.29  FT 

DESIGN  WAVE  PERIOD  =  4.47  SEC 

DESIGN  DEEP  WATER  WAVE  LENGTH  =  102.80  FT 

FETCH  FOR  WIND  SETUP  =  7.50  MILES 

WIND  SETUP  =  0.42  FT 

SIGNIFICANT  WAVE  RUNUP  =  5.46  FT  ABOVE  SWL 

TOTAL  INCREASE  IN  WATER  LEVEL  =  5.89  FT  ABOVE  SWL 

ENTER  RERUN  OR  STOP 
-STOP 


The  following  is  a  list  of  the  permanent  input  data  files  DH7780  and 
DH77801  as  changed  by  the  last  run  of  each  file. 
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100  H7780 
102  MCGEE  BEND 


RE5ERV0IR 


104 

STATION 

A 

106 

15  -1 

10  0 

108 

20.00 

73.00 

2.50 

110 

55.00 

50.00 

112 

1 . 33 

2.70 

2.89 

114 

12.31 

11.74 

6.39 

1.00 

3.22 

5.87 


6.49  12.64  16.95 

5.59  3.46  3.41 


DH77801 


100  •r773E 

1C2  rnr  cT|  TPCfO* 

10/1  '■-p-irr  * 

105  i5-i  i  £  e 

I  £  i  2  0.00  52.32' 

113  55.23  45.00 

112  2.23  2.21 

II  4  7.36  7.  53 


2. 52  1.33 

2.92  3.22  5.  05 

£.11  1.71  1.24 


4.  53 

1.25 
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REF:  ER  1110-1-10  -  ENGINEERING  AND  DESIGN  -  Engineering  and  Computer 
Program  Library  Standards  and  Documentation,  Appendix  C 

PART  III:  FILE  DOCUMENTATION 

1.  REVISION  LOG:  Program  coding  revised  June  78  to  accommodate  the 
CORPS  time-share  features. 

2.  TITLE:  H7780  -  Wave  Runup  and  Wind  Setup  -  Computational  Model 

3.  PROGRAM  SOURCE  LISTINGS:  See  pages  59-68 

4.  NUMERICAL  AND  LOGICAL  ANALYSIS:  Wave  runup  and  wind  setup  are 
solved  by  direct  solution  of  algebraic  equations  except  for  the  transi¬ 
tional  depth  wave  length  which  is  solved  via  an  interactive  technique. 

5 .  SUBROUTINES  NOT  DOCUMENTED  IN  ABSTRACT:  None 

6.  MISCELLANEOUS :  The  program  is  part  of  the  CORPS  computer  system. 
CORPS  is  an  acronym  standing  for  Conversationally  Oriented  Real-Time 
Program-Generating  System.  The  program  is  now  operational  on  the  WES 
G635»  Vicksburg,  MS;  HIS  66/80,  Macon,  GA;  and  Boeing  CDC,  Seattle,  WA. 
The  source  listing  on  page  59  contains  the  first  line  run  command  and 
brief  for  H7780.  This  first  line  run  command  runs  the  binary  H7780B  of 
the  source  listing  on  pages  60-68  (FORTRAN  source  of  H7780)  and  attaches 
the  WESLIB  routines  RERUN,  HACCT,  and  TACHFILE. 
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0001##RUN  WESLIB/CORPS/H7780B,RiWESLIB/RERUN,R;WESLIB/HACCT,R; 
0002*#WESLIB/TACHFILE,R 

0800  62THE  REQUIRED  INPUTS  FOR  THIS  PROGRAM  CONSIST  OF  THE  JOB  TITLE, 
0805  58THE  SITE  ID,  THE  NUMBER  OF  RADI AL S (15  MINIMUM),  THE  RADIAL 
0810  63LENGTHS  IN  MILES,  WHETHER  TO  USE  DESIGN  OVERLAND  WINDSPEED  OR  1 
0815  63MIN  AND  60  MIN  WIND  SPEED,  THE  1MIN  AND  60MIN  WIND  SPEED  OR  THE 
0820  63DESIGN  OVERLAND  WIND  SPEED  IN  MPH,  THE  FACTOR  FOR  A  SIGNIFICANT 
0825  62DESIGN  WAVE  OR  FOR  OTHER  RELATED  DESIGN  WAVE.  WHETHER  TO  PRINT 
0830  62THE  DEEP  WATER  WAVE  CHARACTERISTICS  OR  NOT,  THE  WATER  DEPTH  AT 
0835  61THE  TOE  OF  THE  STRUCTURE  IN  FT,  THE  AVERAGE  WATER  DEPTH  ALONG 
0840  60WIND  FETCH  IN  FT,  THE  COTANGENT  OF  THE  ANGLE  BETWEEN  THE  EM- 
0845  63BANKMENT  SLOPE  AND  THE  HORIZONTAL,  THE  FACTOR  FOR  REDUCING  WAVE 
0850  62RUNUP  ON  SMOOTH  SLOPE  WHERE  THE  SLOPE  IS  SLIGHTLY  ROUGHER  THAN 
0855  6 1SMOOTH ,  THE  SLOPE  COVER! RIPRAP ,  SMOOTH,  OR  OTHER), AND  WHETHER 
0860  55TO  SAVE  THE  TABULATED  OUTPUT  FOR  GRAPHICS  OR  OTHER  USE. 

0865  530UTPUT  INCLUDES  THE  BASIC  DATA  FROM  THE  INPUT  AND  THE 
0870  62COMPUTATIONAL  RESULTS  WHICH  INCLUDE  THE  EFFECTIVE  WAVE  FETCH 
0875  6  0  IN  MILES,  THE  WATER  WAVE  CONDITIONS  USED  FOR  PRESENT  PROBLEM 
0880  6 1 ( TYPE  EMBANKMENT  SLOPE,  DESIGN  WIND  SPEED  IN  MPH,  DESIGN  WIND 
0885  6 1  DURATION  IN  MIN,  DESIGN  WAVE  HEIGHT  IN  FT,  DESIGN  WAVE  PERIOD 
0890  59IN  SEC,  DESIGN  DEEP  WATER  WAVE  LENGTH  IN  FT,  FETCH  FOR  WIND 
C895  62SETUP  IN  MILES,  WIND  SETUP  IN  FT,  SIGNIFICANT  WAVE  RUNUP  IN  FT 

0900  62ABOVE  SWL ,  AND  TOTAL  INCREASE  IN  WATER  LEVEL  IN  FT  ABOVE  SWL ) , 

0905  6 1  AND,  IF  REQUESTED  IN  INPUT,  THE  WIND  DATA  AND  DEEP  WATER  WAVE 

0910  62CHARACTERISTICS(0VER  WATER  WIND  SPEED  IN  MPH,  WIND  DURATION  IN 

0915  63MIN,  SIGNIFICANT  WAVE  HEIGHT  IN  FT,  AND  SIGNIFICANT  WAVE  PERIOD 
0920  6 0 IN  SEC)  AND  WAVE  HEIGHTS  AND  WAVE  RUNUP  FOR  WAVES  OTHER  THAN 
0925  63THE  SIGNIFICANT  WAVECWAVE  EXCEEDANCE  IN  PERCENT,  WAVE  HEIGHT  IN 
0930  22FT ,  WAVE  RUNUP  IN  FT). 

0999#06FINISH 
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00001* #RUN  *=; WESL IB/CORPS/H7780BC  NOGO ) 

09000  CHARACTER  HFILE*5 

09010  COMMON  /MA IN/LQZ , LQX/C7 7805/HFI L E 
09020  HFILE=5HH7780 
09030  LQZ=1;LQX=1 
09040  15000  PRINT  10006 

09050  10006  FORMAT (/"INPUT  H7780  -  WAVE  RUNUP  AND  WIND  SETUP , COMPUTATIO 
09060JNAL  MODEL"//) 

09070  CALL  HACCT (HFILE) 

09090  IFCLQZ.NE.l)  GO  TO  15005 
09100  CALL  H7780($15011) 

09110  GO  TO  15010 
09120  15005  PRINT  9125 

09125  9125  FORMAT ( "CHANGE  SPECIFIC  ITEMS  ON  EXISTING  DATA  SET  OR  NEW  DA 
09127JTA  SET"/"NEW  Y  OR  N") 

09130  CALL  ANSWER($15006, $15007) 

09140  15006  CALL  DET ACH ( 1 , ,  ) ; CAL L  H77803 C $1 50 1 1 ) ; GO  TO  15010 
09150  15007  CALL  H7 780 1 < $  150 1 1 ) 

09160  CALL  H77802($15011) 

09170  15010  LQZ=2 
09180  15011  PRINT,"  " 

09190  CHARACTER  ZZZZZZ*2 

09200  16000  PRINT,  "ENTER  RERUN  OR  STOP" 

09210  READ  16001,  ZZZZZZ 
09220  16001  FORMAT ( A2 ) 

09230  I F( ZZZZZZ . EQ . 2HRE)  GO  TO  15000 
09240  I F ( ZZZZZZ . EQ . 2HST )  GO  TO  20000 
09250  PRINT, "ERROR  ###  RETYPE" 

09260  GO  TO  16000 
09270  20000  STOP;END 
09999  SUBROUTINE  H7780C*) 

10000  COMMON  /MAIN/LQZ,LQX/C77801/KFILE,FILEK(2)/C77802/NSTART,NSAVE,TIT 
10010*LE(15),SITEC15),IMM,ITYPE,IWIND,IDES,ITAB,D1,D,COTTH,FACT,U2,U61,U 
10020*L/C77803/X(50)/C77804/UF(50) ,T(50) , H C 50 ) , TDUR C 50 ) , RP ( 20  ) , HP ( 20  ) ,PP 
10030JEPC20) , IX/C77805/HFILE 
10040  DIMENSION  AL PHA ( 8 ) , FI < 7 ) , WR C 7 ) , GDC  5 ) 

10050  CHARACTER  FI L EK*8 , T ITL E*4 , SI TE*4 , FI L ET*8 , HFI L E*5 

10  060  DATA  FI ,WR,GD,GRAV,SCOSAL ,PI/0 . 5, 1. , 2., 3., 4., 5., 1. , 1.08, 1.13, 1.2 
10070*1,1.26,1.28,1 .3,6.5882, .0161,. 36 92, 2. 2024,. 87 98, 32. 2, 13. 510917 39, 3 
10075* . 14159265/ 

10078  C-22./15. ;RAD=PI/180. 

10080  DO  10090  1=1,8 

10090  10090  ALPHA(I)=(48-I*6)*RAD 

10100  ENTRY  H77803C*) 

10110  NSTART=0;LQZ=1 ;KFILE=1 
10120  PRINT, "PERM  DATA  FILE;Y  OR  N" 

10130  CALL  ANSWER($10200, $10140) 

10140  10140  PRINT, "WE  WILL  HELP  YOU  SETUP  YOUR  DATAFILE.  ENTER  FILE  NAME 
10150*" 

10160  READ.FILEKCKFILE) 
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10170  CALL  H7780I($10190);DS=D1;U1=U2;U60=U61 

10180  GO  TO  10370 

10190  10190  RETURN  1 

10200  10200  PRINT, "DATAFILE  NAME" 

10210  READ,FILEK(KFILE);NSTART=1;LQZ=2 
10220  ENTRY  H7 7 80 2 (x ) ; DS  =  D1 ; U1 =U2 ; U6 0 =U6 1 
10230  10230  CALL  TACHFI L E( HFI L E , $  1 0 1 90  > 

10240  REWIND  KFILE 

10250  READCKFILE, 10260 )  FILET 

10260  10260  FORMAT ( 4X , A8 ) 

10330  10330  READ(KFILE, 10350 , END=10370)  TI TL E , SI T E , IMM , ITYPE , IWIND , IDES , 
103404ITAB. D1,D, COT TH, FACT, U2,U61»(X(J),J=1, IMM );D$=D1 ; U1=U2; U60=U6 1 
10350  10  350  FORMAT (4X, 15 A4/4X, 1 5 A4/4X, 5 1 3/4X , 4 F7 . 2/4X , 2F7 . 2 , 7 ( /4X , 8F7 . 2 ) 
10360*1 

10370  10370  IFCNSTART.EQ.0.AND.LQZ.EQ.2)  GO  TO  10420 
10375  PRINT  10380 

10380  10380  FORMAT ( "CHANGE  ANY  DATA  BEFORE  RUN.Y  OR  N") 

10390  CALL  ANSUERC$10400  ,  $10420  ) 

10400  10400  L  QZ-2 ; CAL  L  H7  78  0  1  <  S 1 0 1 90  ) 

10410  GO  TO  10230 

10420  10420  I F( IWIND . EQ . 1 )  GO  TO  10440 
10430  UL  =U1 ; U6 0  =  0  . 

10440  10440  PRINT  1 0450 , T IT L E , SITE 

10450  10450  FORMAT (///"OUTPUT  H7780  -  WAVE  RUNUP  AND  WIND  SETUP , COMPUTA 
1 04  55*  T ION AL  M0DEL"//15A4/15A4//28X, "BASIC  DATA"//) 

10460  PRINT  10470, D1,D,C0TTH, FACT 

10470  10470  FORMAT! 5X.29HDEPTH  AT  TOE  OF  EMBANKMENT  =  ,F7.2,3H  FT/5X.22H 
1 0480* AVERAGE  WATER  DEPTH  =  ,F7.2,3H  FT/5X, 1 9H EMBANKMENT  SLOPE  =  ,F7.2,4 
10490* HH:1V/5X,35H SMOOTH  SLOPE  MODIFICATION  FACTOR  =  ,F7.2) 

10500  IFCIWIND.EQ.l)  GO  TO  10540 
10510  PRINT  10520, UL 

10520  10520  FORMAT ( 5X, 31 HOB SERVED  OVERLAND  WIND  SPEED  =  ,F7.2,4H  MPH) 

10530  GO  TO  10570 

10540  10540  PRINT  10550, U2.U61 

10550  10550  F0RMAT(5X,24H0NE  MINUTE  WIND  SPEED  =  ,F7.2,4H  MPH/5X , 26HS IXT 
1 0 56 0* Y  MINUTE  WIND  SPEED  =  ,F7.2,4H  MPH) 


10570 

10570 

FE 

=0 . ; 11=1; 

IJK 

*15 

10580 

10580 

J  = 

0 ; IFIN=0; 

SXA 

= 0 . 

10590 

DO  10660 

I=II,IJK 

10600 

I F ( I  FIN . 

EQ.l)  GO 

TO 

10640 

10610 

J  =  J  +  1 

10620 

I  F(  J  . 

EQ. 

8)  I  FIN  =  1 

10630 

GO  TO 

10 

650 

10640 

10640 

J- 

J-l 

10650 

10650 

XA 

=X(I)#COS 

(AL 

PHA( J))**2 

10660 

10660 

SXA=SXA+XA 

10670 

F=SXA/SC 

OSAL 

10680 

IF(  FE 

.GT 

.F)  GO  TO 

10700 

10690 

I  K  =  1 1 

;  fe 

=  F 

10700 

10700 

IFCIJK. EQ. IMM) 

GO  TO  10720 

6i 
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10710  II=II+l;IJK-IJK+l;GO  TO  10580 
10720  10720  PRINT  10730 

10730  10730  F0RMAT(y/8X,58HRADIALS  CORRESPONDING  TO  CRITICAL  EFFECTIVE 
1  0740*  FETCH  IM  MI L ES// 1 
10750  Il-IK+14 

1 Q  76  0  PRINT  10770, CXC I > , I =IK , IL > 

10770  1  0770  FCP,riAT(5X,8F7.2) 

10780  PRINT  10790, FE 

10790  10790  F0RMATC///25X, "COMPUTATIONAL  RESULTS"//5X, 23HEFFECTIVE  WAVE 
10800JFETCH  =  ,F7.2,6H  MILES) 

10810  FES=FE*5280 . ;RMAX=0 . 

10820  DO  10840  1=IK,IL 

10830  IFCRMAX.LT. X(I))  RMAX=XCI) 

1C840  10840  CONTINUE 

10850  IFCFE.GE.5. )  GO  TO  10930 

10860  INC-0 

10870  10870  INC:INC+1 

10880  IFCFICINC) -GE.FE)  GO  TO  10900 

10890  GO  TO  10870 

10900  10900  RG=(FE-FI(INC-1))/(F1CINC)-FI(INC-I)) 

10910  R=RG* CWRC INC) -WR( INC-1 ) )+UR( INC-1) 

10920  GO  TO  10940 
10930  10930  R - 1 . 3 

10940  10940  IFCIWIND. EO. 0)  GO  TO  11260 

10950  U1-RXU1 

10960  U120=.96*R*U60 

10970  U360=  .8S*R*U60 

10980  U60=U60*R 

10990  IX  =  'J1-U360 

11000  U-AINT (U360) 

11010  A60-ALOGC6C. );A360=ALOGC360. ) ; AU60=AL0GC U6 0 ) ; AU36 0 = ALOGC U36 0 ) 
11020  US=U1*C 

11030  11030  GFU-GRAV*FES/US**2;GF1=AL0GCGFU) 

11040  ALN-GD(2)*GF1XX2 
11050  BLN-GDC 31XGF1-GDC4) 

11060  0LN=GDC5)*GF1 

11070  TUF=(US#GD( 1 ) *EXP ( SORT ( A L N-BLN ) +DLN ) >/ ( GRAVX6 0 . ) 

11080  IFCTDF.LE.60. )  GO  TO  11140 
11090  A1 =( AU60-AU360 )/C A60-A360) 

11100  A0-EXPCAU6Q-A1XA60) 

11110  A1-ABSC1 ./Al) 

11120  TDC-(C*A0/US)##A1 

11130  GO  TO  11190 

11140  11140  U30=.3*U1+.7*U60 

11150  IFCTDF. LE. 30 . )  GO  TO  11180 

11160  TDC  =  CU30-US/C)*30 ./CU30-U60  >  +  30 . 

11170  GO  TO  11190 

11180  11180  TDC=CU1-US/C)«29./CU1-U30)+1. 

11190  11190  IFCTDF. GT.TDC)  GO  TO  11210 
11200  US=US+.l;GO  TO  11220 
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11210  11210  US=US-1. 

11220  11220  1 r ( TDF-TDC . LT . .  2 )  GO  TO  11240 

11230  GO  TO  11030 

11240  11240  IUD=US/C 

11250  UD=IUD; GO  TO  11290 

11260  1 i 26  0  IX  =  20 

11270  UD=R*UL ; I UD=UD; UD= IUD 

11230  U-UD-1 0 . 

11290  11290  AU  =  0. 

11300  DO  11450  1=1, IX 
11310  UF( I ) =U+ AU 
11320  US=UF(I)*C 

11330  GrU=GRAV*FES/US*X2;GFl=AL0G(GFU) 

1  1340  GT=TANH( . 077XGFU** -251 
11350  GH=TAHH( . 0125*GFU** .42) 

11360  T(I)=2.4*PI*GT*US/GRAV 
11370  H(I)=.283*US*#2*GH/GRAV 
11380  ALN=GD(2)*GF1**2 
11390  BLN=GD(3)*GF1-GDC4) 

1  1400  DLN  =  GD( 5 ) *GF 1 

11410  TDURC I )=(US*GD(1 ) *EXP ( SQRTC ALN- BLN ) +DLN ) )/ ( GRAV*6 0 . > 

11420  IFCUFC I >  —  UDJ  11440,11430,11440 

11430  11430  TD=T(I);HD=H(I);UD=UF(I); TDURD= TDURC I ) 

11440  11440  AU  =  AU  +  1  . 

11450  11450  CONTINUE 

11460  11460  DI=GRAV*TD*X2/(2.*PI) 

11470  CO  =DL/TD 

11480  IFCD/DL . LE. .5)  GO  TO  11530 
11490  PRINT  1 1 50 0 

11500  11500  F0RMATC///8X.51HDEEP  WATER  WAVE  CONDITIONS  USED  FOR  PRESENT 
1 1 5 1  OiPROF!  L  EM ) 

11520  GO  TO  11850 

11530  11530  IFCD/DL. LT. .04)  GO  TO  11580 
11540  PRINT  11550 

11550  11550  F0RMATC///8X.59HTRANSITI0NAL  WATER  WAVE  CONDITIONS  USED  FOR 
11560$ PRESENT  PROBLEM) 

11570  GO  TO  11610 
11580  11580  PRINT  11590 

11590  11590  F0RMATC///11X,54HSHALL0W  WATER  WAVE  CONDITIONS  USED  FOR  PRES 
1 16008 ENT  PR03LEM) 

11610  11610  US=UD*C 
11620  GFU=GRAV#FES/US  *#2 
11630  GDU  =  GRAWD/US**2 
11640  G2=TANHC .833*GDU*#.375) 

11650  G1=TANHC .53*GDU**.75) 

11660  G3= . 283*G1#TANHC . 0 125*GFU** . 42/G1 > 

11670  HD=US**2*G3/GRAV 

11680  G4=l .2*G2*TANHC . 07 7 *GFU** . 25/G2 ) 

11690  TD=2.*PI*US*G4/GRAV 
11700  IFCD/DL. LT. .04)  GO  TO  11820 


63 


H7780 


11710  Al EN=DL 

11720  GTA:GRAV#TD##2/(2.*PI) 

11730  GTB=2.*PI#D 

11740  11740  BLEN=GTA*TANH(GTB/ALEN) 

11750  DIFF-AIEN-BLEN 
11760  ADI FF  =  AB5( DI FF  ) 

11770  IFCADIFF. LE. 1 .  )  GO  TO  11810 
11780  IF(DIFF)  11800,11800,11790 
11790  11790  ALEN=ALEN-1 . ; GO  TO  11740 
11800  11800  ALEN-ALEN+1 . ;GO  TO  11740 
11810  11810  DL  =  AL£N ; GO  TO  11850 
11820  11820  TD-2.*PI*US*G4/GRAV 
11840  DL=TD#SQRTCGRAV*D) 

11850  11850  FS-2 . * FE 

11860  IF(FS. LE.RMAX)  GO  TO  11880 

11870  FS=RMAX 

11830  11880  S=FS*UD**2/(1400.*D) 

11890  DS-DS+S 

11900  IF(ITYPE)  11910,12050,12310 
11910  11910  IF(C0TTH.LT.5.  )  GO  TO  11960 
11920  PRINT  11930 

11930  11930  F0RMAT(/1X,57HPR0GRAM  CANNOT  COMPUTE  RUNUP  ON  SLOPES  FLATTER 
11940XTMAN  1  ON  5) 

11950  GO  TO  12310 

1  1960  11960  1 r ( COTTH . GT . 2 .  )  GO  TO  12010 
11970  PRINT  11980 

11980  11980  FORMAT (/70HPROGRAM  CANNOT  COMPUTE  RUNUP  ON  A  VERTICAL  OR  NEA 
11990XR  VERTICAL  EMBANKMENT) 

12000  GO  TO  12310 

12010  12010  RS=HD/(.4+SQRT(HD/(GRAV#TD*#2/(2.*PI)))*COTTH) 

12020  PRINT  12030 

12030  12030  F0RMAT(/2X,25H EMBANKMENT  SLOPE  -  RIPRAP//) 

12040  GO  TO  12300 

12050  12050  IF( COTTH- 1 .5)  12060,12210,12060 
12060  12060  IF(C0TTH.GE.2.25)  GO  TO  12110 
12070  FRINT  12080 

12030  12030  FORMAT (/IX, 49HPR0GR AM  CANNOT  COMPUTE  WAVE  RUNUP  FOR  SLOPE  GI 
12090XVEN) 

12100  GO  TO  12310 

12110  12110  IF(COTTH.GT.6. )  GO  TO  12190 
12120  ZA-C1 .58-2. 35/COTTH ) *SQRT ( HD/DS ) 

12130  ZB-.092*COTTH-.26 
12140  ZF-ZA+ZB 

12150  SINB-1  . /SQRT(C0TTH##2+1.  ) 

12160  ZC=5IN3<(5.95/C0TTH+1.5) 

12170  RS-HD#ZC«(  .123*DL/HD)#*ZE 
12180  GO  TO  12230 

12190  12190  RS=.4#TD‘5QRT(GRAV#HD)/ COTTH 
12200  GO  TO  12230 

12210  12210  ZE=.56*SQRT(HD/DS)-.18 
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12220  RS=HD#2.3*( .123XDL/HD)*#ZE 
12230  12230  R$=FACT*RS 
12240  I F( FACT-  1  .  )  12280,12250,12280 
12250  12250  PRINT  12260 

12260  12260  FORMAT ( /2X , 25HEMBANKMENT  SLOPE  -  SMOOTH) 

12270  GO  TO  12300 
12280  12280  PRINT  12290 

12290  12290  FORMAT ( /2X , 34HEMB ANKMENT  SLOPE  -  MODIFIED  SMOOTH) 

12300  12300  RSS  =  S  +  RS ; GO  TO  12340 
12310  12310  PRINT  12320 

12320  12320  FORMAT ( /5X , 51HRUNUP  IS  TO  BE  DETERMINED  FROM  FIG.13.ETL  1110 
123308-2-21 ) 

12340  12340  PRINT  12350 , UD, TDURD , HD, TD , DL , FS , S 

12350  12350  FORMAT ( /5X, 20HDESIGN  WIND  SPEED  =  ,F7.2,4H  MPH/5X, 23HDESIGN 
123608  WIND  DURATION  =  ,F7.2,4H  MIN/5X, 21 HDESI GN  WAVE  HEIGHT  =  ,F7.2,3H 
123708  FT/5X, 21HDESIGN  WAVE  PERIOD  =  ,F7.2,4H  SEC/5X, 3 2H DESIGN  DEEP  WATE 
1 2 38 0 8 R  WAVE  LENGTH  =  ,F7.2,3H  FT/ 5X , 23HFETCH  FOR  WIND  SETUP  =  ,F7.2,6H 
123908MILES/5X, UNWIND  SETUP  =  , F7 . 2 , 3H  FT) 

12400  IF(ITYPE.GT.O)  GO  TO  12440 
12410  PRINT  12420, RS.RSS 

12420  12420  FORMAT < 5X , 25HSIGNI FI  CANT  WAVE  RUNUP  =  ,F7.2,13H  FT  ABOVE  SWL 
124308/5X, 32HT0TAL  INCREASE  IN  WATER  LEVEL  =  ,F7.2,13H  FT  ABOVE  SWL) 
12440  12440  I F( I  TAB . EQ . 0 )  GO  TO  12610 
12450  PRINT  12550 

12550  12550  FORMAT (/// 14X , 45HWIND  DATA  AND  DEEP  WATER  WAVE  CHARACTERISTI 
1 256 08C3//6X, 10 HOVER  WAT ER , 8X , 4HWIND , 2X , 2 C 6X, 1 1HSIGNI FIC ANT )/6X , 1 OHWI ND 
125708SPEED,6X,8HDURATI0N,6X, 11HWAVE  HEIGHT , 6X, 11HWAVE  PERI0D/8X, 5HCMPH 
125808), 11X,5H( MIN), 10X,4H(FT),13X,5H(S EC)//) 

12590  PRINT  12600, ( UF( I ) , TDURd ) , H( I ) , TCI), 1=1, IX) 

12600  12600  FORI1ATCC6X,F7.2,9X,F7.2,8X,F7.2,10X,F7.2)) 

12610  12610  IF(IDES.EQ.Q.OR.ITYPE.GT.Q)  GO  TO  12740 
12620  PRINT  12630 

12630  12630  FORMAT (///IX, 70HWAVE  HEIGHTS  AND  WAVE  RUNUP  FOR  WAVES  OTHER 
126408  THAN  THE  SIGNIFICANT  WAVE/ 3X , 6 5H ( NOTE-  THE  TERM  WAVE  EXCEEDENCF 
126508  REFERS  TO  THE  PERCENT  OF  WAVES  IN/5X,44HA  WAVE  SPECTRUM  THAT  EXCE 
1 266 08 EDS  A  GIVEN  VALUE. )//lX, 2 ( 12X , 4HWAVE ) , 8X, 4HWAVE/ 10X, 10HEXCEEDANCE 
1267088X,6HHEIGHT ,6X,6H  RUNUP/ 1  OX , 1 OHIN  PERCENT , 9X , 4H ( FT ), 8X , 4H ( FT )// , 
12680  DO  12710  1=1,20 
12690  AI=I;P=A1'100.;PPEP(I)=AI 
12700  RPC  I )=RS*SQRT< ALOG( 1  ,/P)/2.  ) 

12710  12710  HPCI)=HD*RPCI)/RS 

12720  PRINT  12730  , ( PP EP ( I ) , HP ( I ) , RP ( I )  , 1  =  1, 20) 

12730  12730  FORMAT ( ( 1  OX , F7 . 2 , 9X , F7 . 2 , 5X , F7 . 2  ) ) 

12740  12740  I F ( NSA VE . NE . 1 . AND . NS A VE . NE . 2 )  GO  TO  12770 
12750  KFI l E  =  2 ;  CALL  H7 78 OH ( NS AVE ) 

12760  K F I L E  =  1 
12770  12770  RETURN 
12780  END 

14000  SUBROUTINE  H7780K*) 

14010  COMMON  /MAIN/LQZ, L QX/C7 780 1/KFI L E, FI L EK C 2 )/C77802/NSTART , NSAVE , TIT 
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14020*LE(15),SITE(15),IMM,  ITYPE,IWIND.IDES,ITAB,D1.D»C0TTH»FACT,U2,U61,U 

14030SL/C7780  3/X( 5  0  VC77S0 5/HFILE 

14040  CHARACTER  F I L  EK*£  , TITLE*4,SIT£#4,HFILE*5 

14050  I F ( L  QZ . EQ . 3 )  GO  TO  14090 

14060  ENTRY  H77801OO 

14070  IF(LQZ.EQ.2.0R.NSTART.EQ.O)  GO  TO  14090 

14080  L  QZ  - 2  ;  REIJI ND  KFILE 

14090  14090  KKK-14 

14100  IF(LQZ.EQ.l)  GO  TO  14140 

14110  14110  CALL  RERUN(KKK, LQX, JKL) 

14120  GO  10(14140,14180,14210,14240,14270,14310,14390,14440,14480,14510, 
14130*1454 0,14530, 14630, 14660, 14730), JKL 

14140  14140  PRINT , "AA-ENTER  JOB  TITLE  <  OR  =  60  CHARACTERS" 

14150  READ  14160, TITLE 
14160  14160  FORMAT ( 15A4 ) 

14170  GO  TO  (14180, 14110), LQZ 

14180  14180  PRINT, "AB-ENTER  SITE  ID  <  OR  =  60  CHARACTERS" 

14190  READ  14160, SITE 
14200  GO  TOU4210, 14110), LQZ 
14210  14210  PRINT  14215 

14215  14215  FORMAT ("AC-ENTER  THE  NUMBER  OF  RADIALS  USED  TO  DETERMINE  EFF 
14217XECTIVE"/" FETCH  DISTANCED  MINIMUM  OF  15  RADIALS  IS  REQUIRED") 

14220  READ , IMM 

14230  GO  TO(14240, 14110), LQZ 
14240  14240  PRINT  14245, IMM 

14245  14245  FORMAT ( "AD-ENTER  THE  ",I2,"  RADIAL  LENGTHS  IN  MILES, SPEARATE 
1 424 7*  BY  COMMAS") 

14250  READ, (X(I), 1=1, IMM) 

14260  GO  T0(14270, 14110), LQZ 

14270  14270  PRINT  14280  _  „ 

14230  14280  FORMAT( "AE-USE  DESIGN  OVERLAND  WIND  SPEED  OR  1  MIN  AND  60  MI 
1 4 2 9 0 * N  WIND  SPEED'V'AS  IN  ETL  1 1 1 0-2-22 1 ; 0 , DESIGN ; 1 , 1  MIN  AND  60  MIN") 
14300  READ, IWIND 

14310  14310  IF(IWIND.NE.l)  GO  TO  14350 

14320  PRINT, "AF-ENTER  1  MIN  AND  60  MIN  WIND  SPEEDS, MPH” 

14330  READ,U2,U61 
14340  GO  TO  14380 

14350  14350  PRINT , "AF-ENTER  DESIGN  OVERLAND  WIND  SPEED, MPH" 

14360  READ, UL 
14370  U2  =  L".  ;U61-0. 

14380  14330  GO  TO ( 1439 0 , 141 1 0 ) , LQZ 
14390  14390  PRINT  14400 

14400  14400  FORMAT ("AG-ENTER  FACTOR  FOR  DESIGN  WAVE ; 0 , DESIGN  WAVE  IS  SIG 
14410*NIFICAN1 "/"WAVE; 1 , DESIGN  WAVE  IS  OTHER  RELATED  WAVE") 

14420  READ, IDES 

14430  GO  TOU  4440, 14110)  ,  LQZ 

14440  14440  PR  I  NT , "AH-ARE  DEEP  WATER  WAVE  CHARACTERISTICS  TO  BE  PRINTED; 
14450*0, NO  PRINT;1, PRINT" 

14460  READ, ITAB 

14470  GO  TOU4480, 14110),  LQZ 
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14480  14480  PRINT, "AI-ENTER  WATER  DEPTH  AT  TOE  OF  STRUCTURE, FT" 

14490  READ , D1 

14500  GO  T0(14510, 14110), LQZ 

14510  14510  PRINT, "AJ-ENTER  AVERAGE  WATER  DEPTH  ALONG  WIND  FETCH, FT" 
14520  READ, D 

14530  GO  TOC14540, 14110), LQZ 

14540  14540  PRINT , "AK-ENTER  COTANGENT  OF  ANGLE  BETWEEN  EMBANKMENT  SLOPE 
14550SAND  HORIZONTAL" 

14560  READ.COTTH 

14570  GO  10(14580, 14110), LQZ 

14580  14580  PRINT14590 

14590  14590  FORMAT ( "A L -ENTER  FACTOR  FOR  REDUCING  WAVE  RUNUP  ON  SMOOTH  SL 
14600SOPE  WHERE  SLOPE  IS"/"SL IGHTLY  ROUGHER  THAN  SMOOTH") 

14610  READ, FACT 

14620  GO  TOU4630, 14110  ),  LQZ 

14630  14630  PRI NT , "AM-ENT ER  SLOPE  COVER  ; -1 , RIPRAP ; 0 , SMOOTH; 1 , OTHER" 

14640  READ, ITYPE 

14650  GO  T0(14660, 14110), LQZ 

14660  14660  PRINT , "AN-SAV E  TABULATED  OUTPUT  FOR  GRAPHICS  OR  OTHER  USE" 
14670  CALL  ANSWER ( $1 468 0 , $  1 47 7 0  ) 

14680  1468C  NSAVE-1 ;KFILE=2 
14690  PRINT, "FILE  NAME" 

14700  READ,FILEK(KFILE) 

14710  CALL  TACHFIIECHFILE, $14760) 

14720  14720  GO  TO  ( 147 30 , 1 4 1 1 0 ) , LQZ 

14730  14730  K FI L E-l 

14740  CALL  H7780W($14/60, $14750) 

14750  14750  NSTART=0 ; RETURN 
14760  14760  RETURN  1 
14770  14770  NSAVE=0;GO  TO  14720 
14780  END 

16000  SUBROUTINE  H7780W(*,x) 

16010  COMMON  /MAIN/ LQZ, LQX/C7780 1/KFIL E , FILEK(2)/C77802/N5TART,NSAVE»TIT 
16020«LE(15),SITE(15),IMM,ITYPE, IWIND , IDES , IT AB , D1 , D, COTTH , FACT , U2 , U6 1 , U 
16030XL/C77803/X(50)/C77805/HFILE 
16040  DIMENSION  LINE(13) 

16050  CHARACTER  FILEK*8,TITLE*4,SILE*4,HFILEX5 
16060  CALL  TACHFILECHFILE, $16240) 

16070  REWIND  KFILE 
16080  LINE! 1 )  =  100 
16090  M8=IMM/8;NLIN=M8+6 
16100  MR=IMM-M8*8 
16110  IF(MR.EQ.O)  GO  TO  16130 
16120  MS=M8+1 ;NLIN=NLIN+1 
16130  16130  DO  16140  I=2,NLIN 
16140  16140  LINE(I)=LINE(I-l)+2 

16150  WRITE! KFILE, 16 180)  L INEC 1 ) , L INEC2 ) , TITLE, LINE( 3 ) , SITE, L INE( 4 ) , IMM, 
161608IT YPE, IWIND, IDES, ITAB,LINE(5),D1,D, COTTH, FACT, LINE(6),U2,U61, (LINE 
16170* ( 1+6 ),( X( J ), J =( 1-1 )*8+l, IX8-(8-MR)*( I/MS)). 1=1, M8) 

16180  16180  FORMAT ( 13, "  H7 78 0" , 2 ( /1 3 , IX , 15 A4 )/ 13 , IX, 51 3/13 , IX, 4F7 . 2/1 3 , 1 
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I6190*X,2F7.2,7(/I3,1X,8F7.2)) 

16200  CALL  DETACH ( KFILE , .  ) 

16210  I F  ( NST  ART . EQ . 1 )  GO  TO  16250 

16220  N5TART=0 

16230  RETURN  2 

16240  16240  RETURN  1 

16250  16250  NSTART=0;  RETURN 

16260  END 

20000  SUBROUTINE  H7780HI NSA VE ) 

20013  COMMON  /C7 780 1/KFI L E , F I L EK( 2 ) /C7 7804/UF ( 50 ) , T ( 50 ) , H ( 50 ) , TDUR( 50 ) , R 

20020*P(20),HP(20),PPEP(20),IX 

20030  CHARACTER  FILEKv8,AFILE*10»FORM#7 

20040  DIMENSION  LINE(50) 

20050  IFCHSAVE.NE.l)  GO  TO  20260 
20060  WRITE(KFILE, 20070)  (IX, 1=1, 4) 

20070  20070  FORMATC’IO  H7780  07  07  EDGE"/"11  ( 07 (/ ) , ( 3X, 7 ( IX, F7 . 2 ) ) ) "/"l 
20 080 *2", 4 (IX, 13) , 3 ( 2X , "20" )/" 1 3" , 4 ( 2X , "1" ) , 3 ( 2X , "5" ) ) 

20090  WRITE(KFILE, 20100) 

20100  20100  FORMAT ( " 14  CURVE  DESIGNATIONS  FOR  H7780  ARE:"/”15  l=OVER  WA 
20110STER  WIND  SPEED" , 5X , ”2=141 ND  DURATION"/"16  3 =SI GNI FICANT  WAVE  HEIGH 
20120ST  4  =  SIGNI FICANT  WAVE  PERIOD"/"17  5=WAVE  EXCEEDANCE" , 1 IX, "6 =WAVE 
20130*  HEIGHT"/"18  7=WAVE  RUNUP"/"19  UNITS  FOR  ABOVE  VARIABLES  ARE:"/"2 
20140*0  MPH  =  1  ITIN-2  SEC=4  FT=3,6,7  PERCENT=5") 

20150  20150  LINE(1)=100 

20160  DO  20170  1=2, IX 

20170  20170  LINE(I)=LINE(I-l)+2 

20172  IF(IX.GE.20)  GO  TO  20180 

20174  DO  20176  I=IX+1,20 

20176  20176  L INE( I ) =L INE( 1-1 )+2 

20180  20180  WRITE(KFILE, 20200)  ( LINE( I ) , UF( I ) , TDUR( I ) , H( I ) , T( I ) ,PPEP( I ) , 
20190*HP(I),RP(I),I=1,20) 

20200  20200  FORMAT ( ( 13 , 7 ( IX , F7 . 2 ) )  ) 

20210  IF(IX.LE.20)  GO  TO  20240 

20220  WRITEOCFILE, 20230)  ( L IHE( I ) , UF( I ) , TDUR( I ) , H ( I ) , T ( I ) , I =21 , IX) 

20230  20230  FORMAT ( ( 1 3 , 4 ( IX , F7 . 2 ) , 4X , "0 . " , 2 ( 6X, "0 . " )  ) ) 

20240  20240  CALL  DET ACH ( KFI L E , , ) ; NS AVE=2 
20250  RETURN 

20260  20260  ENCODE ( AFI LE , 2027 0 )  FILEK(KFIIE) 

20270  20270  FORMAT ("/", A8 ) 

20280  CALL  ATT ACH ( KFI L E , AFI L E , 3 , 0 , , ) ; LH=1 1 
20290  EHCODECFORM, 20300)  LH 
20300  20300  FORMAT ("(", 12 ,"(/))" ) 

20310  READ(KFILE. FORM) 

20320  GO  TO  20150 
20330  END 
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APPENDIX  B 

Design  Wind  Speed  and  Wind  Duration 
Determinations 

The  method  used  in  the  present  computer  program  for  determining 
the  design  wind  speed,  based  on  regional  wind  data,  is  for  all  practical 
purposes  consistent  with  the  approach  presented  in  ETL  1110-2-221.  The 
method  presented  in  the  ETL  requires  plotting  two  wind  velocity-duration 
curves,  one  based  on  regional  wind  and  duration  data,  and  one  based  on 
data  from  wave  forecasting  curves  for  a  given  wave  fetch.  Example  wind 
velocity-duration  curves  are  shown  in  Figure  D-l  of  the  ETL.  The  inter¬ 
section  of  the  two  curves  are  taken  to  yield  the  design  wind  speed  and 
corresponding  wind  duration.  In  the  present  scheme,  equations  are 
written  to  describe  the  wind  velocity-duration  curve  based  on  regional- 
wind  data.  The  empirical  formula  used  for  constructing  wind  velocity- 
duration  curves  in  connection  with  the  wave  forecasting  curves  is 
equated  to  the  equations  describing  the  wind  velocity-duration  based 
on  regional  data.  By  an  iterative  process  this  technique  yields  the 
design  wind  and  corresponding  wind  duration.  The  procedures  and 
formulae  used  for  determining  these  design  values  are  covered  in  the 
following  paragraphs . 

Three  equations  are  used  to  reproduce  the  wind  velocity-duration 
curve  for  the  regional  wind  data.  One  equation  provides  solutions  for 
wind  speeds  with  durations  from  1  to  30  minutes,  one  for  wind  speeds 
with  durations  from  30  minutes  to  60  minutes,  and  one  for  wind  speeds 
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with  durations  from  60  to  360  minutes.  As  discussed  in  paragraph  5-C-l, 
basic  wind  input  consists  of  wind  speeds  corresponding  to  wind  durations 
of  1  and  60  minutes.  The  wind  speeds  for  these  durations  are  defined 
symbolically  as  U  and  UgQ  ,  respectively.  Equations  of  straight 
lines  are  used  to  obtain  solutions  of  wind  duration  as  a  function  of 
wind  speed  for  wind  durations  ranging  from  1  to  30  minutes  and  from  30 


to  60  minutes.  The  30-minute  wind  speed,  U^q  ,  is  estimated  by 


U30  = 


4*)- 


(U1  -  U60} 


(B-l ) 


The  factor  K  is  taken  to  be  0.U0  based  on  the  curve  shown  in  Fig¬ 
ure  D-l  in  the  ETL.  Thus  this  number  is  simply  a  reduction  factor  for 

estimating  the  30-minute  wind  speed.  The  wind  duration  t.  ,  for  the 

dc 

interval  between  1  and  30  minutes  as  a  function  of  the  wind  speed  is 
given  by 


t ,  =  1  + 

dc 


29  (Ux  -  U) 
U1  '  U30 


1  i  *dc  ±  30 


(B-2) 


and  for  the  interval  between  30  and  60  minutes  by 


30  (U  -  U) 

t .  =  30  +  — r: - — -  ;  30  <  t  <  60 

dc  ~  dc  - 


( B— 3 ) 


30 


6o 


where  U  is  the  actual  wind  speed  in  the  specific  wind  duration 
interval. 

In  the  wind  duration  interval  between  60  and  360  minutes  a 
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logarithmic  transformation  is  used  to  describe  wind  duration  as 
tion  of  the  wind.  This  relation  is  given  by 


a  func- 


dc 


exp(ln  UgQ  -  A  In  60) 


U 


1/A 


60 


<  t  < 
—  dc  — 


360 


(B-l 


where  A  is 


U60  -  In  U360 
In  60  -  In  360 


The  symbol  exp(x)  =  ex  and  In  represents  the  natural  logarithm. 

The  design  wind  and  corresponding  wind  duration  can  be  found 
equating  either  equation  B-2,  B-3,  or  B-1*  to  equation  (M  by  varying 
the  wind  speed  U  in  an  appropriate  range.  This  is  done  by  an 


iterative  process  in  the  computer  routine. 


